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ABSTRACT 

We present a study exploring a systematic effect on the brightness of type la supernovae using numerical 
models that assume the single-degenerate paradigm. Our investigation varied the central density of the progen- 
itor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production 
and distribution of radioactive 56 Ni, which powers the light curve. We performed a suite of two-dimensional 
simulations with randomized initial conditions, allowing us to characterize the statistical trends that we present. 
The simulations indicate that production of Fe-group material is statistically independent of progenitor central 
density, but the mass of stable Fe-group isotopes is tightly correlated with central density, with a decrease in 
the production of 56 Ni at higher central densities. These results imply progenitors with higher central densities 
produce dimmer events. We provide details of the post-explosion distribution of 56 Ni in the models, including 
the lack of a consistent centrally-located deficit of 56 Ni, which may be compared to observed remnants. By 
performing a self-consistent extrapolation of our model yields and considering the main-sequence lifetime of 
the progenitor star and the elapsed time between the formation of the white dwarf and the onset of accretion, 
we develop a brightness-age relation that improves our prediction of the expected trend for single degenerates 
and we compare this relation with observations. 

Subject headings: hydrodynamics — nuclear reactions, nucleosynthesis, abundances — supernovae: general 
— white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNela; singular SNIa) are bright, tran- 
sient astronomical events identified by a peak-light spectrum 
showing no evidence of hydrogen but absorption lines of 
singly-ionized silicon (Minkowski 1941; Filippenko 1997). 
These events follow from explosive thermonuclear burning of 
degenerate stellar material composed principally of C and O, 
which synthesizes ^0.6M Q of radioactive 56 Ni. The decay of 
this 56 Ni powers the light curve (Truran et al. 1967; Colgate 
& McKee 1969; Arnett 1982; Pinto & Eastman 2000). 

The progenitor systems of these explosions remain the sub- 
ject of considerable debate and active research. Observa- 
tions, however, indicate these events largely form a homo- 
geneous class. Phillips (1993) identified a relationship be- 
tween the maximum B-band magnitude of an event and its 
rate of decline. This "brighter equals broader" relationship 
has been extended to additional bands with templates from 
nearby events, allowing these events to be calibrated as an 
extension of the astronomical distance ladder (see Jha et al. 
2007 for a description). This property, along with the bright- 
ness of SNela, which makes them visible over great distances, 
enables the use of SNela to probe the structure and expan- 
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sion history of the universe, allowing studies of various cos- 
mological models' parameters (Riess et al. 1998; Perlmutter 
et al. 1999; Albrecht et al. 2006; Kirshner 2010), with recent 
work constraining cosmological parameters to within a few 
percent (Riess et al. 2011; Sullivan et al. 2011). Recent ob- 
servational studies of SNela have begun to correct for cor- 
relations of the brightness of a SNIa with properties of the 
host galaxy (Conley et al. 2011). Many SNIa observations 
are restricted to broadband photometry, so knowledge of host 
galaxy properties is correlated. The inability to deconvolve 
these properties from each other is among the larger sources 
of uncertainty in cosmological constraints from SNela, so ad- 
vancing the understanding of how brightness correlates with 
host galaxy properties may contribute significantly to reduc- 
ing the uncertainties of cosmological parameters. 

The brightness, and therefore "broadness", of a SNIa is de- 
termined principally by the amount of 56 Ni synthesized dur- 
ing the explosion. Observations report that SNela appear to 
have an intrinsic scatter of a few tenths of a magnitude after 
calibration, forcing a minimum uncertainty in any distances 
measured by using SNela as standardizable candles (Jacoby 
et al. 1992; Kirshner 2010). An important goal of theoreti- 
cal research into SNela, from the standpoint of cosmology, is 
to understand the sources of scatter and to identify potential 
systematic biases by studying the effects of various proper- 
ties on the mechanism and nucleosynthetic yield of the SNIa. 
The surrounding stellar population, the metallicity and mass 
of the progenitor, the thermodynamic state of the progenitor, 
the cooling and accretion history of the progenitor, and other 
parameters are known to affect the lightcurves of SNela; the 
role of these "secondary" parameters is the subject of con- 
siderable study (e.g., Ropke et al. 2006; Htiflich et al. 2010). 
Additionally, many of these effects may be interconnected in 
complex ways (Domfnguez et al. 2001; Lesaffre et al. 2006; 
Townsley et al. 2009). 
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Observational campaigns are gathering information about 
SNela at an unprecedented rate. Scannapieco & Bildsten 
(2005) and Mannucci et al. (2006) showed that the delay 
time (elapsed time between star formation and the supernova 
event) data are best fit by a bimodal delay time distribution 
(DTD) with a prompt component that tracks less than 1 Gyr 
after star formation and a tardy component that occurs sev- 
eral Gyr later. Gallagher et al. (2008) demonstrate a correla- 
tion between brighter SNela and shorter delay times, which 
they state is consistent with the bimodality described by Man- 
nucci et al., but also with a continuous relation. Howell et al. 
(2009), Neill et al. (2009) and Brandt et al. (2010) also find 
such a correlation between the delay time and brightness of a 
SNIa. While the degeneracy of age and metallicity in obser- 
vations could obscure these correlations, Howell et al. (2009) 
note that the scatter in brightness of this observed relation is 
unlikely to be explained by the effect of metallicity. 

For this theoretical study, we adopt the model known as 
the single-degenerate paradigm. This model assumes that a 
SNIa is the result of a thermonuclear disruption of a white 
dwarf (WD) in a mass-transferring binary system with either a 
main-sequence or red-giant companion star (see Branch et al. 
1995; Filippenko 1997; Hillebrandt & Niemeyer 2000; Livio 
2000; Ropke 2006; Li et al. 201 1; Nugent et al. 201 1; Bloom 
et al. 2012, and references therein). Recent observational evi- 
dence, however, suggests other progenitors such as the merg- 
ing of two white dwarfs may explain many events (Scalzo 
et al. 2010; Yuan et al. 2010). In the single-degenerate sce- 
nario, the WD is formed when the primary star goes through 
a giant phase and expels a planetary nebula. Once the pri- 
mary becomes a WD, it is initially not in contact with the 
companion star, and it slowly cools as thermal energy is ra- 
diated away. Once the companion star evolves and fills its 
Roche lobe, mass-transfer begins to carry low-mass elements 
from the envelope of the companion to the surface of the WD. 
If the accretion rate exceeds ~40~ 7 Mg yr , the H-rich mate- 
rial can steadily burn (Nomoto et al. 2007) and the WD gains 
mass, which heats and compresses the WD, driving up both 
the temperature and density in the core. Once the temperature 
rises enough for carbon burning to begin, the core of the WD 
begins to convect; this is known as the "simmering" phase. 
This simmering phase lasts on order of 10 3 yr, and ends when 
a flame is ignited, which occurs approximately when the eddy 
turnover time becomes shorter than the local nuclear runaway 
time. Our initial models attempt to parameterize the WD at 
the end of the simmering phase, just at the beginning of the 
thermonuclear deflagration, which will in turn cause an ex- 
plosion that will disrupt the entire WD in a SNIa. 

The explosion mechanism we use (within the single- 
degenerate paradigm) is that of a deflagration to detonation 
transition (DDT). After ignition, the flame propagates as a 
subsonic deflagration for a while and then transitions to a su- 
personic detonation that rapidly consumes the star (Blinnikov 
& Khokhlov 1986; Woosley 1990; Khokhlov 1991; Hoflich 
et al. 1995; Hoflich & Khokhlov 1996; Khokhlov et al. 1997; 
Niemeyer & Woosley 1997; Hoeflich et al. 1998; Niemeyer 
1999). We describe the details of our implementation of this 
explosion mechanism below. 

In the single-degenerate paradigm, a longer delay time can 
be explained by a longer elapsed time between the formation 
of the WD and the onset of accretion. During this period, 
the WD is in isolation and cools, hence the moniker the "WD 
cooling time" (t coo i)- Following the cooling time is a period 
of accretion, during which the WD is compressed and heats, 



approaching the conditions for ignition of the thermonuclear 
runaway. The decrease in temperature during the cooling 
time, which is determined by r coo i, influences the density 
structure of the WD just prior to ignition, with a longer r coo i 
resulting in a higher central density when the core reaches the 
ignition temperature (Lesaffre et al. 2006). Thus, a correlation 
between central density and the brightness of an event would 
suggest a correlation between delay time and the brightness 
of an event. 

In this manuscript, we expand on our earlier investiga- 
tion on the effect of r coo i on the brightness of the explosion. 
In Krueger et al. (2010) we reported that as the central density 
of the progenitor WD increases, the production of radioactive 
56 Ni decreases due to increased neutronization rates, produc- 
ing a dimmer event. Using the results of Lesaffre et al. (2006), 
we related the WD central density to r coo i and were able to 
compare our results to the observations of Neill et al. (2009). 
Here we present additional details of our models; a statisti- 
cal analysis of the results including the assessment of intrin- 
sic scatter; the distribution of Fe-group elements within the 
remnant; and a potentially-observable effect to demonstrate 
the connection between age, progenitor central density, and 
brightness. We also revised our previously-reported trend in 
brightness with age to account for the main sequence evolu- 
tion of the WD progenitor. 

In Section 2 we discuss the methodology of our suite of 
simulations, followed by details of the code we used to per- 
form our simulations in Section 3. We present the results of 
our simulations in Section 4, and discuss how these results 
compare with previous studies in Section 5. Section 6 con- 
tains a brief summary and final conclusions. 

2. METHODOLOGY 

As described above, our explosion models assume the DDT 
mechanism in which a flame ignited in the core propagates as 
a subsonic deflagration and then transitions to a detonation. 
We simulate an explosion from ignition through the detona- 
tion phase until burning effectively ceases. Our models are by 
necessity incomplete, however, in that we do not produce light 
curves and spectra with which we could compare to actual ob- 
servations. Instead, we rely on the mass of 56 Ni synthesized 
in our models and compare this result from our simulations to 
56 Ni masses inferred from observations (Howell et al. 2009). 

For the study, we adopted the theoretical framework first 
applied in Townsley et al. (2009) for a statistical study of a 
suite of simulations performed with randomized initial condi- 
tions. We constructed a set of five progenitor models with dif- 
ferent central densities (p c ,o), and therefore slightly different 
masses, and from these performed suites of two-dimensional 
simulations. For each progenitor model we applied thirty sets 
of initial conditions consisting of randomized perturbations 
on an initially-burned region. We refer to each of these thirty 
as a realization, each of which is seeded by a random number 
used to generate a unique power spectrum of spherical har- 
monics (see Appendix A for details of the realizations). The 
spectra are used as initial perturbations to a spherical "match 
head" in the center of the progenitor star. Each progenitor 
WD had the same seed values applied, resulting in the same 
thirty perturbations. This procedure allows us to characterize 
the intrinsic scatter in the models and check for systematic 
biases in the realizations across different progenitors, such as 
how the morphology of the initial conditions may influence 
the final result. 

Complete details of our models and simulations follow, but 
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we preface the description by mentioning that one limitation 
of our study is the use of two-dimensional models and the 
parameterization of inherently three-dimensional phenomena. 
In particular, two-dimensional models lack any meaningful 
consideration of turbulence and its effect on the flame because 
the turbulence found in the interior of a simmering white 
dwarf and its interaction with the flame are inherently three- 
dimensional. The problem is compounded by the fact that this 
interaction occurs partially on unresolvable scales, necessitat- 
ing use of sub-grid-scale models (see Schmidt et al. 2006a, 
for an example). 

There are two critical ways that omitting consideration of 
turbulence influences this work. The first is in the calcula- 
tion of flame speeds during the deflagration phase. The burn- 
ing model (described below) relies on an input flame speed 
to propagate a model flame during the deflagration. Turbu- 
lence will interact with this flame, stretching it and thereby 
boosting the burning rate (see Schmidt et al. 2006a,b, and ref- 
erences therein). Our present models boost the input flame 
speed from tabulated laminar values to compensate for buoy- 
ancy effects (Townsley et al. 2009) but do not include the 
effect of turbulence-flame interactions for reasons discussed 
above. Effectively, we assume the increase in flame surface 
(and hence the burning rate) is dominated by stretching due to 
buoyancy rather than turbulence. 

The second critical way that omitting consideration of tur- 
bulence influences this work is in the criteria for the DDT (see 
Seitenzahl et al. 201 1, and references therein). For this study, 
we parameterized the DDT criterion as a threshold density, 
Pddt, with a detonation initiated when the top of a rising 
plume of burned material reaches this density. As we de- 
scribe below, this threshold density determines the duration 
of the deflagration and thus strongly influences the outcome 
of an explosion. Our a priori choice for this parameter led 
to higher than expected yields of radioactive 56 Ni, necessitat- 
ing a rescaling of our results for comparison to observations. 
Confirmation of the trends we present from extrapolated re- 
sults awaits a future study with a more consistent treatment of 
these issues. 

2.1. Initial White Dwarf Models 

As demonstrated by Lesaffre et al. (2006), for a given zero- 
age main-sequence mass, the properties of a progenitor WD 
such as central temperature and density at ignition of the de- 
flagration can be constructed as functions of t coo i. We take the 
leading-order effect from varying t coo i; that is, vary p c fi while 
holding all other parameters constant. This choice allows us 
to disentangle the effects of p c o from the effects of other pa- 
rameters. For our p c o, we chose 1 -5 x 10 9 g cm -3 in steps of 
1 x 10 9 g cm" 3 . We then constructed a series of five parame- 
terized WD progenitor models in hydrostatic equilibrium. 

Figure 1 presents the profiles of the progenitor WDs in the 
p-T plane. The core of each WD is isentropic due to convec- 
tion, the intermediate ("envelope") region is isothermal due to 
the high conductivity of degenerate matter, and the outer ("at- 
mosphere") region has a power-law temperature dependence 
that was chosen to mimic a radiative atmosphere. For a model 
to explode, the central temperature must be in the range where 
the carbon burning begins a runaway, which is approximately 
7-8 x 10 8 K. Varying the isothermal envelope temperature 
would have a relatively insignificant effect on the mass, as the 
envelope contains only a small fraction of the mass; primarily 
the mass is set by p c i q. Thus we have chosen a central tem- 



TABLE 1: Central densities, masses, and radii of progenitor 
WDs. 



Pcfi (g cm "') 


Mm < M ©) 


M core (M ) 


R(km) 


1 x 10 9 


1.345 


1.180 


2500 


2 x 10 9 


1.368 


1.162 


2076 


3 x 10 9 


1.379 


1.144 


1852 


4x 10 9 


1.385 


1.131 


1716 


5 x 10 9 


1.389 


1.121 


1604 



TABLE 2: Composition of the progenitor WDs. 



isotope 


mass fraction 


core envelope 


12 C 


40% 50% 


16 Q 


57% 48% 


22 Ne 


3% 2% 



perature at the low end of the carbon ignition range and allow 
the total mass to vary as a function of p c o. Table 1 shows the 
total mass (M tot ) and mass of the isentropic core (M core ) for 
each progenitor. 

The core of each WD has a lower C/O ratio than the en- 
velope; this is primarily due to the composition of different 
regions of the star at the end of the asymptotic giant branch 
(AGB) phase and the subsequent mixing of these regions, with 
additional contributions from the consumption of C during 
the simmering phase (Straniero et al. 2003; Piro & Bildsten 
2008; Chamulak et al. 2008; Piro & Chang 2008). Our pa- 
rameterized models for this study assume a fixed carbon abun- 
dance. As we plan to also study the dependence on central 
carbon abundance, we have maintained a clear separation be- 
tween the central-density and the carbon-abundance studies 
by varying only central density and not the core carbon frac- 
tion. The composition discontinuity between the core and the 
envelope causes a temperature discontinuity (Piro & Chang 
2008), shown by the short vertical line segments between the 
core and the envelope in Figure 1 . The composition is listed 
in Table 2. We use 22 Ne as a placeholder to represent the 
neutron-rich isotopes present in a SNIa. The abundance of 
22 Ne is calibrated to achieve the electron-to-baryon ratio of 
the material present in a SNIa, but sedimentation effects are 
not included. See Jackson et al. (2010), specifically Section 2, 
for a more detailed discussion. 

2.2. Ensemble of Simulations 

The thirty unique realizations for each of our five progenitor 
models, consisting of a set of perturbations on the initially- 
burned region, allowed us to perform a suite of 150 two- 
dimensional, axisymmetric simulations of SNela. A simu- 
lation begins with a region of burned material at the center of 
the star with the perturbation from sphericity given as spher- 
ical harmonics with a prescribed range ofl2<^<16 and 
random amplitudes; the amplitudes for each realization are 
shown in Appendix A. Each realization has a unique seed 
(for the random number generator), allowing each realization 
to be applied to the five progenitor models. Figure 2a shows 
two example perturbations that span the space of the random 
perturbations. Realization 21, on the left, is representative 
of the "spikiest" initial conditions; i.e., the greatest deviation 
from the mean radius. Realization 10, on the right, is repre- 
sentative of the "smoothest" initial conditions; i.e., the small- 
est deviation from the mean radius. 

3. THE SIMULATION CODE 
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Fig. 1. — : The structure of our progenitor WDs in p-T space. 
Below ~ 10 4 g cm" 3 is the power-law atmosphere. In the cen- 
ter of the figure is the isothermal envelope, with the same tem- 
perature for all progenitors. In the upper right region of the 
figure is the adiabatic core; all 5 progenitors are isentropic in 
this region, but the value of the entropy varies between pro- 
genitors. 

The simulations were performed using a customized 
version of the FLASH code 9 , an Eulerian adaptive- 
mesh compressible hydrodynamics code developed by the 
ASC/Alliances Center for Astrophysical Thermonuclear 
Flashes at the University of Chicago (Fryxell et al. 2000; 
Calder et al. 2002). The equation of state (EoS) that we use is 
the fully-ionized electron-ion plasma EoS (Timmes & Swesty 
2000; Fryxell et al. 2000). The version of FLASH used here is 
the same as that used in Jackson et al. (2010). Customizations 
perform two main functions: First to implement the energy 
release due to explosive carbon-oxygen fusion, in both defla- 
gration and detonation propagation modes, as well as a provi- 
sion for a transition from deflagration to detonation. Second, 
criteria for mesh refinement that capture the important physics 
with a suitable degree of efficiency. Details of various com- 
ponents are given by Townsley et al. (2007), Townsley et al. 

(2009) , and Jackson et al. (2010), with important additional 
and supporting information on implementation of some pieces 
of physics in Calder et al. (2007) and Seitenzahl et al. (2009b). 

As mentioned in the introductory material, the simulations 
implement a multi-dimensional version of the deflagration- 
to-detonation transition model, in which the flame is born as 
a subsonic deflagration and later transitions to a supersonic 
detonation (Blinnikov & Khokhlov 1986; Woosley 1990; 
Khokhlov 1991; Hoflich et al. 1995; Hoflich & Khokhlov 
1996; Khokhlov et al. 1997; Niemeyer & Woosley 1997; Hoe- 
flich et al. 1998; Niemeyer 1999). The initial conditions pre- 
scribe the deflagration at the start of a simulation. As the de- 
flagration proceeds, the flame is subject to fluid instabilities 
(Figure 2b), and when the top of a rising bubble reaches the 
threshold density, it is assumed to transition into a supersonic 
detonation (Figures 2c and 2d; see Section 3.1.2 for details 
of the transition). See Townsley et al. (2009), Maeda et al. 

(2010) , Jackson et al. (2010), Ropke et al. (201 1), Seitenzahl 
et al. (2011), and references therein for examples of recent 
work assuming this explosion mechanism. 

3.1. Burning Model 
9 available from http://flash.uchicago.edu 



Since the carbon-oxygen fusion occurring in SNela pro- 
ceeds to nuclear statistical equilibrium (NSE), in which all 
reactions among all nuclides are, to good approximation, fast 
relative to the hydrodynamic timescales, the nuclear process- 
ing necessarily involves a large number of reactions and nu- 
clides. This is especially true when calculating e~ capture on 
Fe-group elements (IGEs), as the overall effective e~ capture 
rate is a combination of contributions from captures on a wide 
variety of nuclides with comparable individual rates. As ex- 
plored in our earlier work (Calder et al. 2007; Townsley et al. 
2007; see also Khokhlov 2000), it is possible to abstract the 
burning process from hundreds of nuclides to just a few fluid 
state variables with appropriately chosen reaction dynamics 
and energetics computed using a large set of nuclides. This 
enables tremendous gains in computational efficiency, mak- 
ing 3-d simulations and extensive 2-d studies such as this fea- 
sible. 

Our burning model consists of three reaction progress vari- 
ables that describe conversion between four states. The first 
state is the initial, unburned mixture of 12 C, 16 0, and 22 Ne 
(representing general neutron excess, as described in section 
2.1), which we call fuel. The second state is the result of 
the 12 C fusing to roughly Si-group elements, which we call 
ash. The third state is the result of the remaining 16 burning 
to Si-group elements, which we call nuclear statistical quasi- 
equilibrium (NSQE) material. The fourth and final state is the 
result of NSQE material relaxing to nuclear statistical equi- 
librium (NSE), containing mostly IGEs. The three reaction 
progress variables are 

4>f a Carbon consumption, fuel to ash 

<t>aq Oxygen consumption, ash to NSQE 

4> qn Conversion of Si-group to Fe-group, NSQE to NSE. 

Each variable evolves from (unburned) to 1 (fully burned). 
We also require that <fif a > (f) aq > cf> q „ to enforce the time- 
ordering of the four states. A given cell will have mass frac- 
tions \ -<f) f a of fuel, 4>f a - 4> aq of ash, (j) aq - <\> qn of NSQE ma- 
terial, and <f) qn of NSE material. 

While a large portion of the energy release occurs in the 
consumption of 12 C, the final NSE state is not, during the ex- 
plosion, energetically inert. The NSE state is one in which the 
distribution of nuclides in the fluid, and therefore the average 
nuclear binding energy, is determined by the fluid state. That 
is, it participates in the EoS of the fluid, releasing or absorbing 
energy as the pressure and density of the fluid change. As a re- 
sult, a significant portion of the energy release for some fluid 
elements occurs well after the "fully burned" NSE state has 
been reached. The reaction kinetics used are given in Jackson 
et al. (2010). We track several material properties in the fully 
burned material, including the electron-to-baryon ratio (Y e ), 
ion-to-baryon ratio (Y-, on ), and average nuclear binding energy 
perbaryon (q). 

Our treatment of matter in NSE allows for the effects of 
weak reactions, specifically electron captures, which serve to 
deleptonize the material. Weak processes (e.g. electron cap- 
ture) are included in the calculation of the energy input rate, 
as are neutrino losses, which are calculated by convolving 
the NSE distribution with the weak interaction cross sections. 
Both the NSE state and the electron capture rates were cal- 
culated with a set of 443 nuclides (Seitenzahl et al. 2009b). 
Weak cross sections were taken from Fuller et al. (1985), Oda 
et al. (1994), and Langanke & Martinez-Pinedo (2001), with 
newer rates superseding earlier ones. 
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Fig. 2. — : Illustrations of the phases of the DDT model of a SNIa; each panel displays two of our realizations: 21 (left) and 10 
(right). Panel 2a shows the initial flame surfaces; realization 21 has the greatest deviation from the mean radius, and realization 
10 has the least deviation from the mean radius. Panels 2b - 2d are snapshots from our simulations, with p c o = 3 x 10 9 g cm" 3 ; 
Panel 2b shows the early deflagration phase, Panel 2c shows the first DDT events, and Panel 2d shows the later detonation 
phase. The colors represent the four stages of the burning model discussed in Section 3.1: cream represents unburned fuel, 
gold represents ash from carbon burning, red represents material in NSQE, and black represents material in NSE. Only in the 
detonation stage do the three burning processes separate out spatially; they are co-located during the deflagration. Our simulations 
extend to ~6.5 x 10 8 km, but these images show only the inner regions; the spatial scale varies between panels, with only 
panels 2b and 2c having the same spatial extent in order to illustrate the expansion that occurs during the deflagration phase. 
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The treatment of electron capture is critical to the dynamics 
for three reasons. First, NSE is a dynamic equilibrium and 
the composition of material in NSE evolves as the thermody- 
namic state evolves and/or Y e changes. Binding energy can be 
released if the equilibrium evolves toward more tightly bound 
nuclei, which changes the local temperature. Second, the re- 
duction in Y e lowers the Fermi energy, reducing the primary 
pressure support of this highly degenerate material and hav- 
ing an impact on the buoyancy of the neutronized material. 
Finally, neutrinos are emitted (since the star is transparent to 
them) so that some energy is lost from the system. 

By using the progress variables defined above and the local 
Y e , we can derive a local estimate for the abundance of 56 Ni. 
As was done in Townsley et al. (2009), we estimate the 56 Ni 
abundance by assuming that the first IGE material made as 
neutronization occurs is equal parts by mass 54 Fe and 58 Ni. 
The local mass fraction of 56 Ni is estimated by 



Ye-(l-$ q n)Y e .f 



Xs6 Ni = max 



0.5-0.48212' 



(1) 



where Y e j is the electron fraction in the unburned fuel and 
0.48212 is the electron fraction of material that is equal parts 
by mass 54 Fe and 58 Ni. This should be a modestly accurate 
estimate of the 56 Ni production because the dynamics of the 
progress variable <fi q „ have been calibrated to reproduce, in 
hydrodynamics, the production of IGEs during incomplete 
silicon burning seen in direct calculations (Zeldovich, Von 
Neumann, Dtiring; ZND; e.g. Khokhlov 1989) of steady-state 
detonations. When realization 2 from Jackson et al. (2010) 
with pddt = 10 7 1 g cm" 3 is post-processed with a nuclear net- 
work, the 56 Ni yield from Equation (1) is within 2% of that 
determined by the post-processing. More detailed study of 
the accuracy of computed yields under various conditions is 
the subject of separate, ongoing work. 

3.1.1. Deflagration 

Even at our highest resolution (4 km), the flame front is 
unresolved. To handle this we use an artificial, resolved reac- 
tion front that is governed by the advection-diffusion-reaction 
(ADR) equation (Khokhlov 1995; Vladimirova et al. 2006), 
with special features to ensure that the front is stable and 
acoustically quiet (Townsley et al. 2007). Our ADR front is 
chosen to be resolved over about 4 computational cells in or- 
der to obtain acceptably low acoustic noise as it propagates 
across the grid and releases energy. This creates an extended 
"partially burned" region that requires some specialized treat- 
ment. In such regions, particularly at high density, we say that 
the material is well-separated into unburned and fully-burned 
material divided by a thin flame. However the spatial resolu- 
tion cannot capture this and the average over a mixture of fully 
burned and unburned results a "partially-burned" state. Thus 
we have to make estimates of the correct thermodynamic state 
of the two cases (unburned and fully burned) mixed together 
within the region. 

Additionally, we enforce a minimum flame speed in order 
to prevent the flame from being torn apart by Rayleigh-Taylor- 
induced turbulence. The minimum flame speed is 



S min = 0.5 y/Agm A, 



(2) 



where A is the Atwood number, g is the local acceleration of 
gravity, A is the width of the grid cell, and m is an adjustable 



parameter, set to 0.04 for these simulations. The flame speed 
is set by 

S = max(S mi n,Siam)- (3) 

The Atwood number and the laminar flame speed, S\a m , are 
both functions of the local, unburned density estimate and the 
composition. The Atwood number varies by less than 0.01% 
due to the amount of 22 Ne present so it is tabulated for a rep- 
resentative, constant 22 Ne fraction. 

3.1.2. DDT 

At present, the physical mechanism by which a DDT in 
degenerate supernova material occurs is an area of current 
research (see Rtipke 2007; Seitenzahl et al. 2009a; Woosley 
et al. 2009; Schmidt et al. 2010; Poludnenko et al. 2011, 
and references therein). Simulations of supernovae involv- 
ing a DDT assume it occurs via the Zeldovich-gradient mech- 
anism (Khokhlov et al. 1997, but see also Niemeyer 1999), 
in which a gradient in reactivity leads to a series of explo- 
sions that are in phase with the velocity of a steadily propa- 
gating detonation wave. Many authors suggest that when the 
flame reaches a state of distributed burning, which is when 
turbulence on scales at or below the laminar flame width are 
fast enough to dominate transport processes (see, e.g., Pope 
1987), fuel and ash are mixed and the temperature of the fuel 
is raised and "prepared" in such a way to produce the re- 
quired reactivity gradient. A requirement for distributed burn- 
ing is that the ratio of turbulent intensity to the laminar flame 
speed must exceed some unknown threshold, which is still 
actively researched (Niemeyer & Woosley 1997; Khokhlov 
et al. 1997; Golombek & Niemeyer 2005; Ropke & Hille- 
brandt 2005; Aspden et al. 2008, 2010; Poludnenko & Oran 
2011a,b). Entrance into the distributed burning regime does 
not guarantee such a reactivity gradient to form. Woosley 
(2007) and Woosley et al. (2009) studied incorporating more 
stringent requirements for these conditions to be met. 

In the context of supernova models, the ratio of turbulent in- 
tensity to laminar flame speed changes most rapidly due to the 
change in laminar flame properties, which are strongly depen- 
dent on fuel density. Therefore, DDT is assumed to occur at a 
range of densities that vary somewhat but generally lie in the 
range of 10 6 7 to 10 7 7 g cm" 3 (Khokhlov et al. 1997; Lisewski 
et al. 2000; Woosley 2007; Ropke & Niemeyer 2007; Bravo 
& Garcfa-Senz 2008; Maeda et al. 2010; Jackson et al. 2010). 
We choose to ignite detonations where the flame reaches a 
specific density, Pddt= 10 g cm" 3 , which puts us some- 
where in the middle of that range, consistent with a number 
of studies on this subject. We note that while parameterizing 
the DDT criteria by only density omits effects such as back- 
ground turbulent intensity and density gradients, our choice 
of one threshold is intended to keep such variables constant 
so that we can isolate and investigate the effect of varying the 
central density. 

The choice of Pddt increases or decreases the duration 
of the deflagration phase, which increases or decreases the 
amount of expansion prior to the detonation and hence the 
yield of IGEs. A companion study of the effects of varying 
Pddt as a proxy for metallicity under the same statistical en- 
semble we use here indicated a slight over-production of 56 Ni 
at this choice of Pddt (Jackson et al. 2010), a result borne 
out by this study. The DDT transition is implemented by 
burning small regions ahead of rising plumes. When a plume 
reaches pddt, a circular region with a radius of 12 km is se- 
lected 32 km radially outward from the point where the rising 
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plume reaches the transition density. The reaction progress 
variables in this region are then instantly increased to a fully- 
burned state. This method conserves energy, as the detonation 
is initiated by the sudden release of energy from the conver- 
sion of fuel to NSE, not by an unphysical addition of extra 
energy. Each plume is limited to ignite no more than 2-3 
detonations, with a minimum separation distance of 200 km 
imposed between detonation points. The full details of this 
algorithm are given in Section 3.2 of Jackson et al. (2010). 

Previous studies using this DDT mechanism found that 
DDT points with a 12 km radius successfully generated det- 
onations in all simulations performed. However, our simula- 
tions showed that this size DDT point is not as robust for this 
study. One of our simulations (p c p = 1.0 x 10 9 g cm" 3 , real- 
ization 8) deviated from the behavior of the other 149 simula- 
tions and inspection showed that the first plume to reach Pudt 
ignited several DDT points that did not propagate as detona- 
tions; the first detonation to actually propagate ignited when 
the second plume reached /?ddt significantly after (and be- 
neath) the first plume to reach pddt- The failure of the det- 
onation of the first plume to reach /?ddt led to significantly 
more expansion of the WD prior to the subsequent detona- 
tion, which led to some or all of the material burning at a 
lower density than it should have, and the corresponding sus- 
pect yields. Thus we removed this simulation from our suite 
and performed the analysis on the remaining 149 points. Even 
if this "failed DDT" simulation is included, the results pre- 
sented in this paper change by no more than a few percent. 

3.2. Mesh Refinement 

Refinement is based on gradients in p and <pf a , subject to 
the limits imposed below. We define three types of regions, 
subject to different refinement limits: 

1. fluff (f): regions with p < p^g 

2. star (*): non-energy-generating stellar material, p > 

Pfluff 

3. energy generation (eg): 

regions with e nuc > e eg or ^ fa > (f)^ 

where </>f a is the reaction progress variable from the 
ADR equation, and e eg and </>f a , eg are parameters equal to 
10 18 erg g -1 s" 1 and 0.2 s" 1 respectively. These limits are 
chosen so that all actively propagating flames or detonation 
fronts are at the highest resolution. We establish a mini- 
mum cell size for refinement of each type of region such that 
Af > A* > A eg . We use A eg = 4 km, A* = 16 km and A/ 
to be as large as allowable. FLASH only allows adjacent sub- 
domains of the mesh to be of fixed size (we use 16x16 cells) 
and to differ by a single refinement level (a factor of 2 in 
resolution). These were found to be the lowest resolutions 
which gave converged results in 1-d simulations (Townsley 
et al. 2009). 

Fluff is the low-density area outside of the star; we choose 
Pfluff = 10 3 g cm -3 . The FLASH code cannot properly handle 
empty (zero-density) regions, so these regions are set to some 
low, but non-zero, density so that they will not affect the dy- 
namics of the star. To avoid rapidly cycling the refinement- 
derefinement of a region, a small amount of hysteresis is in- 
troduced near the limits for refinement changes, so that, for 
example, a refinement of grid resolution is not immediately 
derefined as a result of slight changes due to the necessary 
interpolation. 



4. RESULTS 

The results presented here build on the initial results pre- 
sented in Krueger et al. (2010) and extend the analysis be- 
yond what was shown there. The explosions occur in two 
main phases, the deflagration phase (from ignition until the 
first DDT occurs) and the subsequent detonation phase. The 
principal difference is that during the deflagration phase, the 
star has time to react to the energy release. Accordingly, the 
evolution is naturally divided by the time at which the first 
DDT event occurs, which we define as ?ddt; this duration also 
describes the time spent in the deflagration-dominated phase 
of the SNIa evolution. We also define ?ige as the time when the 
production of IGEs ceases. This is the time at which burning 
ceases, and by this time the NSE state is no longer evolving 
due to freezeout; thus energy release has effectively ceased 
by fiGE- However, our models have not yet entered into free 
expansion by this time. For the purpose of our simulations, 
fiGE also measures the duration of the entire SNIa event. We 
found that the duration of the deflagration phase (?ddt) de- 
creases with increasing density and is less sensitive to density 
as the density increases, while the duration of the detonation 
phase (equal to ?ige — ^ddt) is very nearly constant for all sim- 
ulations, with a mean of 0.476 s and a standard deviation of 
0.065 s. 

The table in Appendix C presents the masses of 56 Ni and 
IGEs at foDT and Tige- Recall from Table 1 that the total mass 
and the mass of the convective core both change with den- 
sity, but the variations are only of a few percent. Thus, we 
believe that effects from variations in the total or convective 
core masses are negligible. The results and trends we describe 
follow from variations in the central density and the variations 
in initial conditions from realization to realization. Also, as 
discussed in Section 3 of Krueger et al. (2010), the choice of 
the DDT transition density in our simulations led to an over- 
production of 56 Ni. Essentially, our models are systematically 
too bright, but we believe that our trends are valid. Jackson 
et al. (2010) investigated the role of DDT density in our mod- 
els and found that the production of 56 Ni is very sensitive to 
the choice of DDT density. This choice determines the dura- 
tion of the deflagration phase, which determines the amount 
of expansion and, accordingly, the density profile of the star 
during the detonation and the yield. Future studies will be 
better calibrated based on these results. 

4.1. Evolution 

Figure 3 shows the evolution of the gravitational binding 
energy and the mass of IGEs for the five simulations per- 
formed using realization 5, each with a different p c $. During 
the deflagration-dominated phase, higher-p e o progenitors ex- 
perience faster burning, thus expanding the star faster due to 
the faster energy release. This effect can be seen by the rapid 
drop of the binding energy for the higher-density simulations 
around 0.5 s in the upper panel. The transition to a detonation 
in higher- p C Q progenitors occurs sooner. Once the first DDT 
event occurs, the production of IGEs proceeds much faster 
than in the deflagration phase, as may be seen by the sudden 
increase of the slope in the lower panel after the detonation 
occurs. 

During the detonation phase the drop in binding energy is 
similar for all densities. During this same period the curves 
of the IGE mass stay ordered, with the differences between 
them following principally from the difference in ?ddt- The 
nonlinear morphological dependencies come into play as the 
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Fig. 3. — : Sample curves showing the evolution of the gravita- 
tional binding energy and the mass of IGEs over the duration 
of a simulation. These curves all use realization 5. Curves 
are colored by p c ,o: 1.0 x 10 9 g cm" 3 (red), 2.0 x 10 9 g cm" 3 
(green), 3.0 x 10 9 g cm" 3 (blue), 4.0 x 10 9 g cm" 3 (magenta), 
5.0 x 10 9 g cm" 3 (cyan). The times ?ddt and ?ige are marked 
by black circles and black diamonds respectively. 

detonation slows and then stops; the total IGE yield plateaus 
to a constant value, but that value is not (for a single real- 
ization) correlated with the central density of the progenitor. 
The leveling off of the mass of IGEs, due to the cessation of 
burning, is apparent in the lower panel. The time fjQE was 
calculated for each simulation by finding the point at which 
the IGE mass changes by less than 0.01% over the preceding 
0.01 s. 

4.2. Statistics 

We showed in Krueger et al. (2010) (see especially Figure 
2) that a single initial morphology is, in general, insufficient 
to capture trends in SNela due to the nonlinearities involved 
in the explosion process. This observation invites the ques- 
tion of how many initial morphologies are necessary to obtain 
statistically-meaningful trends from the simulations. Figure 4 
presents the standard deviation of 56 Ni yield as a function 
of the number of realizations, with the realizations added in 
the order presented in Appendix A. Shown are the standard 
deviations for each of the five central densities of the study. 
The obvious evolution of the standard deviation until approx- 
imately 15 realizations are included in the average shows that 
a statistically meaningful average requires approximately 15 
realizations. From this result, we conclude that our our sam- 
ple of 30 realizations per central density is sufficient to fully 
characterize the statistical trends we present. The choice of 
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Fig. 4. — : Plot of the standard deviation of 56 Ni yield as 
a function of the number of realizations, added in the order 
listed in Appendix A, for each of the five central densities of 
the study. The standard deviation of the 56 Ni mass converges 
around 15 total realizations. 

30 realizations allows us to go somewhat beyond character- 
izing the variation in the sample, allowing us to reduce our 
uncertainty on the sample mean. 

This result, the need for an average over an ensemble of 
simulations for determining statistically-meaningful trends, is 
critical to the analysis of our study. Because the results dis- 
cussed in this paper come from a statistical analysis of our 
simulations, the results presented here generally cannot be ap- 
plied predictively to individual SNela, but only apply statis- 
tically to large sets of SNela. Exceptions to this rule will be 
noted explicitly. 

4.3. Initial Morphology Correlations 

We explored whether properties of the initial flame mor- 
phology correlate with the final result. We find that "spikier" 
initial flame surfaces, that is initial morphologies that appear 
to have a small number of large amplitude perturbations, give 
rise to faster plume growth and less expansion, and, therefore, 
a higher yield of IGEs and 56 Ni. Several quantities were tested 
as measures of the morphology of "spikiness", and each gave 
the same qualitative result. For example, Figure 5 plots the 
final yield of IGEs as a function of the maximum radius of the 
initial flame surface and standard deviation of the initial flame 
radius. 

As illustrative cases, consider Figure 2. The right pan- 
els shows realization 10, p c ,n = 3 x 10 9 g cm" 3 , which is a 
very smooth initial flame surface; this same realization is out- 
lined by a black square in Figure 5. The smooth configura- 
tion results in several rising plumes approaching pddt, and 
these multiple plumes burn a larger fraction of the star and 
release more energy than a single plume would during its rise 
to Pddt- The net effect is an increase in both turn and in 
the amount of expansion at ?ddt, leaving less mass at densi- 
ties high enough to burn to IGEs during the detonation phase. 
Contrast this with the left panels of Figure 2, showing real- 
ization 21, p c .n = 3 x 10 9 g cm" 3 , which is a very spiky initial 
flame surface; this same realization is outlined by a black cir- 
cle in Figure 5. For this case, a single plume is dominant over 
all other features of the initial flame surface and rapidly ac- 
celerates towards pom with little competition. This gives a 
short deflagration phase, and burns a lower fraction of the star 
prior to the detonation (see especially Figure 2c). These two 
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effects result in less expansion at fDDT, leaving more mass at 
a density high enough to burn to IGEs during the detonation, 
resulting in a greater yield of IGEs and 56 Ni. 

In the center of each panel of Figure 5, the trend is not 
so clear. Realization 24 (outlined by black diamonds) is a 
good example of the lack of a clear trend in this intermediate 
range. The initial flame surface for realization 24 has multi- 
ple large spikes, so that it has a larger than average maximum 
radius and radial standard deviation, and would be considered 
a spiky case. However, most of these spikes are of compa- 
rable maximum radius, and none of these features develops 
into a dominant plume, as is typical of the spiky cases; in- 
stead, the plumes grow together and the behavior is like that 
of the smooth cases. Thus realization 24 leads to high values 
of ?ddt and low yields of IGEs and 56 Ni relative to simula- 
tions of similar spikiness. Thus the intermediate range does 
not show a strong trend, while the extreme ranges (multiple 
competing plumes starting from a smooth initial surface, or a 
single dominant plume with no significant competition) show 
the trend more clearly. 

4.4. Yields 

Figure 6 shows the masses of IGEs and 56 Ni produced in 
the 149 simulations that were analyzed, along with the 56 Ni- 
to-IGE mass ratio. The figures plot yield vs. central density 
and also show average yields and standard deviations at each 
central density, with best-fit trend lines. The mass of IGEs is 
consistent with a flat line; i.e., it is independent of p c $. How- 
ever, the mass of 56 Ni decreases with increasing p c o- The 
significant scatter in the two masses is readily apparent; the 
mean standard deviations for the IGEs and 56 Ni masses are 
0.108 M Q and 0.105 M Q , respectively. As discussed in Maz- 
zali & Podsiadlowski (2006) and Woosley et al. (2007), as- 
suming a constant IGE mass and varying the 56 Ni mass pro- 
duces SNela that lie approximately along the observed width- 
luminosity relationship, while the width of the relationship 
allows the IGE mass to vary somewhat (c.f. Figures 15 and 20 
of Woosley et al. 2007). We note that the masses of IGEs and 
in particular the masses of 56 Ni from our simulations are on 
average higher than accepted results for masses synthesized 
in actual SNIa events (see Woosley et al. 2007, and references 
therein). 

The 56 Ni-to-IGE mass ratio decreases with increasing p c o, 
as would be expected from a constant IGE mass and a de- 
creasing mass of 56 Ni. However, unlike the constant stan- 
dard deviations of these two masses, the standard deviation of 
the 56 Ni-to-IGE mass ratio increases with p C Q. The variation 
in the 56 Ni-to-IGE ratio is dominated by variation related to 
neutronization: because neutronization can exaggerate differ- 
ences that arise in the hydrodynamics, we find that the stan- 
dard deviation (variation) increases with the total mass of sta- 
ble IGEs synthesized 

Additionally, we found that the yield from burning during 
the deflagration phase is substantially different from the yield 
during the detonation phase. Figure 7 shows the yields of 
stable IGEs and 56 Ni during the deflagration and detonation 
phases. The trend of increasing stable IGE yield and decreas- 
ing 56 Ni yield with increasing central density is most obvi- 
ous in the deflagration phase yield. We interpret this result 
as following from the fact that the detonation phase involves 
burning at densities that are typically lower than those of the 
deflagration phase due to expansion of the star during the de- 
flagration phase. The neutronization rate increases with den- 
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Fig. 5. — : Plots of the mass of IGEs as functions of the maxi- 
mum initial flame radius and the standard deviation of the ini- 
tial flame radius, both being treated as measures of how spiky 
the initial flame surface is, with simulations colored by p Cj n. 
The black square marks realization 10 and the black circle 
marks realization 21, both with p c x) = 3 x 10 9 g cm" 3 . These 
two simulations are shown in Figures 2a - 2d, and are exam- 
ples of extreme cases: very smooth (realization 10) and very 
spiky (realization 21). Black diamonds mark realization 24, 
an example which runs counter to the general trend, and a 
good illustration of the ambiguity of the intermediate regions 
of these plots. 

sity, so it shows up most strongly in the deflagration phase 
before significant expansion occurs. 

4.5. Distribution o/ 56 Ni 

The principal result from a simulation is the mass of 56 Ni 
synthesized in the explosion, which directly sets the bright- 
ness of an event. The synthesized 56 Ni masses are listed in 
Appendix C. The next question to be answered in the analy- 
sis concerns the distribution of 56 Ni in the remnant. Figure 8 
presents radially-averaged profiles of 56 Ni for three of the five 
central densities (p Ci o = 1 x 10 9 , 3 x 10 9 , and 5 x 10 9 ). Also 
shown are representative results from realizations 4, 12, 13, 
and 29 for comparison. The radial profile of 56 Ni mass frac- 
tion varies between simulations, but certain details are consis- 
tent across the ensemble. The inner 0.8 - 1.0 M Q has a high 
56 Ni mass fraction; we refer to this region as the "plateau", 
although there can be significant variation within this region. 
Outside of the plateau, there is a smooth decrease to a 56 Ni 
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Fig. 6. — : Final yields (at Tige) of our simulations plotted 
against p c $- The black lines are the best-fit trend lines, with 
the averages and standard deviations marked by the circles 
and the vertical error bars. 

mass fraction at or near zero at the surface of the WD; we re- 
fer to this as the "decline" region. It appears that this decline 
region may not be fully relaxed to the final profile at time 
fjOE, but may experience some steepening prior to entering 
the free-expansion phase. Lower-/?,, rj simulations have less 
neutronization and therefore generate more Ni. For these 
simulations the 56 Ni mass fraction in the plateau tends to clus- 
ter near the maximum value (~ 0.9), with some deviations 
down as far as ^0.6. The higher-p e o simulations have more 
neutronization and therefore generate less 56 Ni, so the typical 
56 Ni mass fraction in the plateau shows more variation and 
there exist larger deviations from the mean (down to ^0.2). 

The representative realizations presented in Figure 8 
demonstrate the deviation from mean behavior. In partic- 
ular, the representative realizations show "typical" behavior 
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Fig. 7. — : Plots of the masses of stable IGEs (red pluses), 
56 Ni (blue asterisks), and total IGEs (green crosses) produced 
during the deflagration phase (top panel) and the detonation 
phase (bottom panel). The symbols for each mass have been 
horizontally shifted slightly so the symbols do not overlap. 

(the 56 Ni profile nearly matches the mean profile), over- and 
under-luminous models, and 56 Ni holes (patches within the 
plateau with significantly less 56 Ni than the surrounding re- 
gions). Previous studies have reported the presence of a 56 Ni 
hole in the inner region (Htiflich et al. 2010). The green and 
red curves in the bottom panel of Figure 8 (realization 4) il- 
lustrate such a 56 Ni deficit, although the deficit is offset from 
the center of the WD in the case of the red curve. Most of our 
simulations do not show evidence of this 56 Ni hole in the inner 
region. Only a few simulations have such a feature, with 56 Ni 
holes being more common in simulations with a higher p C Q. 
Turbulent mixing caused by the burning processes breaks the 
symmetries that give rise to a consistent central 56 Ni hole in 
1-d simulations. 

Figure 9 shows the distributions of 56 Ni, stable IGEs, and 
non-IGE material for the four sample cases shown in Fig- 
ure 8c (pep = 5 x 10 9 ). For our results, the mass of stable IGEs 
is approximated by any IGE material that is not 56 Ni. The fig- 
ures present the inner 6000 km of the domain, and for all cases 
the bulk of IGEs appears in the plotted region. The distribu- 
tions indicate that 56 Ni holes can be caused by incompletely- 
burned regions (for example, on the axis of Figure 9c just be- 
low the equator) or by neutronized regions (for example, ad- 
joining the previously-mentioned incompletely -burned region 
of Figure 9c, or the band at a radius of ~ 2000 km in Fig- 
ure 9a). The plume rise, velocity fields, and neutronization 
may be asymmetric, which is seen especially in Figures 9a 
and 9d. The degree of asymmetry observed in our models 
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Fig. 8. — : Radially-averaged profiles of 56 Ni mass fraction 
at three central densities (1 x 10 9 , 3 x 10 9 , 5 x 10 9 g cm" 3 in 
Panels 8a, 8b, and 8c respectively). The black curves show 
the mean profile of all simulations at a given p c o, with the 
grey band showing the standard deviation. The four colored 
curves represent four different realizations: 4 (red), 12 (blue), 
13 (green), and 29 (brown). 

suggests that there may be noticeable line-of-sight effects for 
SNela. 

4.6. Distinguishing Age Among SNela of Equal Brightness 

Can we extract from our results insight into the age of 
a progenitor WD given the brightness of the SNIa? From 



our results, one might conclude that the mass of stable IGEs 
(M s tabie) should increase with increasing p Cj o, and that the re- 
lation should have a large scatter (akin to the scatter in the 
56 Ni and IGE relations). Figure 10 illustrates the relationship 
between p c q and M sta bi e , along with the best-fit trend and the 
scatter (the shaded region shows two standard deviations in 
each direction around the best fit curve through the data). The 
data show an increasing trend as expected, with the best fit 
being 

( P c > \ C!< A /Stable V , A m /Stable \ „-.„ 

(4) 

Unlike the 56 Ni and IGE relations, however, there is a very 
small scatter in this relation for M sta bi e , with a standard devi- 
ation of only 0.167 x 10 9 g cm" 3 . This result is the tightest 
relation in our data, to the point where this relation can be 
meaningfully applied to a single event. This tight relationship 
is unlike all of our other relations, which only apply to the 
statistics of large ensembles. 

Given observations of multiple SNela of the same bright- 
ness and a reliable measure of the mass of stable IGEs from 
the observations, our models predict that we can use Equation 
(4) to determine the relative ages of the progenitors. The re- 
lation provides the central density from observed masses of 
stable IGEs. Lesaffre et al. (2006) presents relations between 
the cooling time of the progenitor and the central density of 
the progenitor at the ignition of the thermonuclear runaway. 
Applying the relations of Lesaffre et al. (2006) to the central 
densities, our result allows determination of the cooling times 
the progenitor WDs experienced. Assuming the progenitors 
had the same main sequence mass, we thus obtain a measure 
of the relative ages of the progenitors. Because this result 
is derived from varying the central density of the progenitor, 
we are implicitly assuming that such a variation of p f .n is the 
dominant effect on the mass of stable IGEs. Future work in 
three dimensions will consider central density variation com- 
bined with other effects that may be related to age. 

Our models, specifically the distribution of IGEs in the ex- 
panding remnant, offer insight into what would be required 
for observing the mass of stable IGEs. While our models ex- 
hibit a mild degree of asymmetry, we find systematic behavior 
of the distribution of heavy nuclei for events of a given bright- 
ness (or given mass of 56 Ni). Figure 1 1 shows profiles of 56 Ni 
and IGE mass fractions in three sets. Each set has a simulation 
at each of the five central densities, but the simulations within 
a set are chosen to have approximately the same integrated 
56 Ni mass. Accordingly, the simulations within a set may not 
be from the same realization. These profiles are generated at 
?ige, when the burning has essentially ceased, but not yet into 
the free-expansion phase. 

The plateau (central region) exhibits a slight systematic be- 
havior in Figure 1 1 . Considering both the curves of 56 Ni and 
IGEs, we note a tendency of a less-well-defined plateau in the 
lowest- 56 Ni-mass set (top panel). The start of the decline re- 
gion occurs at a higher enclosed mass for simulations with a 
higher mass of 56 Ni; the plateau extends farther out. Also, 
there appears to be a mild trend in the plateau region of simu- 
lations within a set: the higher p C Q simulations (magenta and 
cyan curves) tend to show a wider plateau in the IGE mass 
fraction. In addition, we note that within these results, some 
simulations exhibit the 56 Ni hole while some do not. The hole 
may be observed in the drastic decreases of some curves at 
the lowest enclosed masses. 
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Fig. 9. — : Plots of the inner 6,000 km of four simulations at fiGE, showing mass fractions of Ni (red), stable IGEs (blue), and 
non-IGE material (white). These are the same simulations plotted in Figure 8c, with p c .o = 5 x 10 9 g cm" 3 ; from left to right, the 
panels show a non-central 56 Ni hole (realization 4), an overluminous case (realization 12), an underluminous case with a central 
56 Ni hole (realization 13), and a typical case (realization 29). The off-center 56 Ni deficit in realization 4 is due to a band of 
neutronized, stable isotopes caused by multiple plumes that have reached a common radius. The central 56 Ni hole in realization 
13 is partially due to a neutronized region around the equator, and partially due to an incompletely -burned region near the axis, 
below the equator. 

cations of our trends, particularly the decreasing proportion 
of Ni with increasing central density. The principle result 
from one of our simulations is the mass of 56 Ni, and there 
are well-established relationships between the brightness of 
an observed thermonuclear supernova and the mass of 56 Ni 
synthesized during the explosion (see Woosley et al. 2007, 
and references therein). Thus our trends may have implica- 
tions for the brightness of events. In this section, we explore 
these connections by comparing our trends in 56 Ni mass to 
56 Ni masses inferred from observations. 

In Krueger et al. (2010) (see especially Figure 4 of Krueger 
et al. (2010), an updated version of which is shown in Fig- 
ure 14 of this paper), we showed that our results agree with 
the general trend of dimmer SNela from older stellar popula- 
tions by comparing with observations from Neill et al. (2009). 
We improve upon that result with two additions: recalibra- 
tion of the 56 Ni mass produced by our models and clarifying 
the definition of "age" by including a range of main sequence 
lifetimes. Additionally, this section discusses how our results 
compare to other studies considering the effect of p c q on 56 Ni. 

5.1. Recalibration of 56 Ni Yield 

As discussed above, on average our suite of simulations ex- 
hibited an overproduction of IGEs and 56 Ni when compared 
to masses estimated from observations of remnants (Woosley 
et al. 2007). We estimate a correction to our 56 Ni masses, 
the details of which are given in Appendix B but which we 
summarize here, that allows us to more directly compare our 
results to 56 Ni masses inferred from observations. 

The recalibration is based on relationships between model 
parameters and explosion yields found in previous studies. 
Jackson et al. (2010) found that higher values of pddt led to 
increased yields due to less expansion of the star during the 
deflagration phase. Townsley et al. (2009) found a correspon- 
dence between the mass at densities above 2 x 10 7 g cm" 3 at 
?ddt and the mass of IGEs synthesized in the yield (Figure 5). 
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Fig. 10. — : Plot of p c o vs. M sta bi e . The shaded region displays 
two standard deviations around the model (Equation (4)). 

We also may consider the distribution of the synthesized 
heavy elements within the models. Figure 12 shows the dis- 
tribution of 56 Ni vs. stable IGEs for 3 of the simulations from 
Fi gure 11a. The densities of the three simulations (p c fl = 
1 x 10 9 , 3 x 10 9 , and 5 x 10 9 g cm" 3 ) span the range of cen- 
tral density. The trend for more stable IGEs at higher p c $ can 
clearly be seen by the more extensive regions colored in blue, 
indicating higher mass fractions of stable IGEs. These 2-d 
plots also show that as p c o increases, the 56 Ni and stable IGEs 
are more sharply segregated; there are more mixed regions (in 
shades of purple) in the lowest-p f ,o simulation (Figure 12a). 

5. DISCUSSION 

The results presented above in Section 4 follow directly 
from our simulations and are as rigorous as can be within 
the limitations of our models (as described in Section 2). We 
may extend our results some in order to investigate the impli- 
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Fig. 11. — : Comparisons of radial profiles of 56 Ni and IGE 
mass fractions for models with nearly the same total mass of 
56 Ni, but different values of p Cj n. The upper set of profiles has 
a mass of 0.689 ± 0.003 Mq, the central set of profiles has a 
mass of 0.763 ± 0.003 Mq, and the lower set of profiles has a 
mass of 0.889 ±0.005 Mq. Lines are colored by p c0 : 1.0 x 
10 9 g cm" 3 (red), 2.0 x 10 9 g cm" 3 (green), 3.0 x 10 9 g cm" 3 
(blue), 4.0 x 10 9 g cm" 3 (magenta), 5.0 x 10 9 g cm" 3 (cyan). 



Although the correspondence found by Townsley et al. (2009) 
was not perfect, these two relationships form the basis for our 
re-scaling. We note that the rescaling assumes that even for 
lower values of pddt, the mass of IGEs is independent of p C Q\ 
this assumption is discussed in greater detail in Appendix B 
based on the expansion characteristics present in our simu- 
lations. However, this is a critical assumption in our anal- 
ysis and is based on a somewhat restricted parameter space 
and therefore is subject to future investigation to confirm it. 
Recent work by Seitenzahl et al. (201 1), who performed full 
three-dimensional simulations covering the range of observed 
56 Ni masses, suggests that this assumption does not hold. See 
section Section 5.3. 

Performing this recalibration lowers the 56 Ni masses (see 
Figure 13), giving yields that are consistent with yields in- 
ferred from observations. A similar trend of decreasing 56 Ni 
mass with increasing central density is also found in these 
extrapolated results. Future models will extend to three di- 
mensions and may include more physically-motivated initial 
conditions or more detailed burning models that capture in- 
teractions between the turbulent velocity field and the flame 
structure. Such changes will require that we re-evaluate our 
choice of model parameters, such as pddt- 

5.2. Comparison to Observations 

One of the principal metrics that we extract from each sim- 
ulation is the mass of 56 Ni produced, which is directly related 
to the brightness of an event. Using the method outlined in 
Howell et al. (2009), we converted stretch reported in obser- 
vational results to 56 Ni masses to compare to the 56 Ni masses 
(derived in Section 5.1) obtained from our simulations. We 
then combined our central density values with the results of 
Lesaffre et al. (2006), which correlate the central density at 
the time of the ignition of the flame front to the cooling time 
of the progenitor WD, allowing us to express our results as 
ages. The results of Lesaffre et al. (2006) suggest that a WD 
with a central density of 1 x 10 9 g cm" 3 will not ignite with- 
out further accretion, so for this comparison we neglect our 
simulations with that value of p c p. Our "age" was previously 
defined as the delay time (r coo i), while the observational re- 
sults used the delay time, which includes the main-sequence 
lifetime (tms)- We can improve our comparison by applying 
a shift to our data to account for tms- We assume that our 
progenitors differ solely in their cooling times in order to sep- 
arate out other effects; therefore we take tms to be constant 
across our results. Our best-fit line would now become 

M 5 6 Ni = alog 10 (r coo i + T M s) + ^, (5) 

where a and (3 are fitting parameters. Because the addition 
of tms is inside the logarithm, this corresponds not only to 
a shift but also an increasing slope. We have selected two 
estimated limiting values for tms (0.05 and 1.0 Gyr, corre- 
sponding to main-sequence masses of approximately 8.0 and 
1.5 M©respectively; see Hansen et al. 2004) and included 
them with the original tms = result in Figure 14. Adding 
in a tms consistent with our C/O progenitors brings our re- 
sults into better agreement with the two right-most points of 
Neill et al. (2009). 

As seen in Figure 14, our theoretical results are not in com- 
plete agreement with the observed data. Observationally, the 
age-brightness correlation may flatten at young ages, while 
our data do not, resulting in our data being overluminous rel- 
ative to young SNela. This study isolated the effects of cen- 
tral density and related that to age assuming a constant main- 
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Fig. 12. — : Comparison of 56 Ni and stable IGE distributions at ?ige for simulations with the same total 56 Ni mass. These 
images are generated from the same simulations used in Figure 11a. From left to right, the central densities are 1.0, 3.0, and 
5.0 x 10 9 g cm" 3 . Red represents the mass fraction of 56 Ni, blue represents stable IGEs, and white represents non-IGEs. Axes 
are in kilometers. As p C Q increases (with constant 56 Ni), more IGEs are produced and less material is mixed (seen by the purple 
regions). 



sequence mass, but there are other effects that may be cor- 
related. Examples of such potentially correlated effects are: 
main-sequence mass and it's correlation with central density, 
metallicity of the progenitor, core 12 C fraction prior to igni- 
tion of the deflagration, sedimentation, and others. Inclusion 
of such effects may modify the results presented here and are 
the subjects of future work. 

5.3. Comparison to Other Theoretical Efforts 

Recent theoretical work by other research groups also ad- 
dresses the role of central density in the single-degenerate pic- 
ture of SNela. Fisher et al. (2010) note that in the gravitation- 
ally confined detonation model (GCD; Plewa et al. 2004), a 
higher central density leads to increased energy release dur- 
ing the deflagration phase, which leads to increased expansion 
of the star and an increase in the production of intermediate- 
mass elements and the corresponding decrease in IGEs. These 
results are consistent with our findings concerning the choice 
of Pddt influencing the production of IGEs discussed above 
in Section 3.1.2. These results are also consistent with some 
of our realizations, but Fisher et al. (2010) do not consider the 
role of neutronization, preventing comparison to our princi- 
pal result that the relative proportion of 56 Ni decreases with 
higher central density, producing a dimmer event. 

Seitenzahl et al. (201 1) performed a study of the DDT sce- 
nario with three-dimensional simulations and a description of 
the flame energetics that accounts for neutronization. Seiten- 
zahl et al. (201 1) similarly find that the 56 Ni fraction of IGEs 
decreases in WD models with higher central densities, but 
they also find that the mass of IGEs increases at higher den- 
sities. The net effect is a roughly constant mass of 56 Ni syn- 
thesized during the explosion, implying that central density 
influences the brightness of an event only as a secondary pa- 
rameter. Seitenzahl et al. (201 1) comment on differences be- 
tween their results and our earlier results (Krueger et al. 2010) 
and, likewise, we offer discussion here. 

First, we note that our results are more similar than they 
might at first appear. Our trend of decreasing 56 Ni follows 
from the increased rate of neutronization at high densities of 



a roughly constant mass of IGEs synthesized during the ex- 
plosion. In a sense, our trend is consistent with that of Seiten- 
zahl et al. (2011). For a given mass of IGEs, the fraction 
that is 56 Ni is lower in WD models with a high central den- 
sity. The difference between our results and Seitenzahl et al. 
(201 1) follows from the the increase in the mass of IGEs syn- 
thesized in their models. If our explosions produced more 
IGEs for higher central density progenitors, then we may well 
see a roughly constant mass of 56 Ni. This observation should 
be readily apparent by considering Figure 6. If the flat slope 
of the IGE yield were instead rising enough, the decreasing 
slope of the 56 Ni yield would be instead rising. Accordingly, 
the key to the differences in our results is understanding the 
reason for the difference in the production of IGEs. 

While there are many differences between the methodol- 
ogy of Seitenzahl et al. (2011) and ours, including (as they 
mention) differences in the flame model (level set vs. ADR), 
differences in the energy release scheme, and structure of the 
computational mesh, perhaps the most substantive difference 
is the dimensionality of the simulations. Three-dimensional 
simulations meaningfully describe turbulent flow, which en- 
ables use of turbulence-flame interaction (TFI) models. Tur- 
bulence in two-dimensional simulations, however, has very 
different properties; particularly, it has an inverse-cascade of 
energy from small to large scales (see Chapter 10 of Davidson 
2004, and references therein). Because of the large scope of 
this study, only two-dimensional simulations were possible. 
As described above, our models use a flame speed that com- 
pensates for buoyancy effects to prevent the flame from being 
torn apart by Rayleigh-Taylor-induced turbulence. Because 
of the issue of turbulence in two-dimensional simulations, we 
chose not to include models for unresolved turbulence and the 
turbulence-flame interaction in our models for this study. 

In their three-dimensional simulations, Seitenzahl et al. 
(2011) similarly account for buoyancy effects, but also im- 
plement a method for accounting for turbulent energy on un- 
resolved scales and the corresponding enhancement to the 
flame speed. Originally proposed by Niemeyer & Hillebrandt 
(1995) and developed in detail by Schmidt et al. (2006a,b), 
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Fig. 13. — : Final yields of our simulations after recalibration 
using Equations (B2) and (B3), the parameter values given 
in Table B2, and the time shifts given in Table B3. The black 
lines are the best-fit trend lines, with the averages and standard 
deviations marked by the circles and the vertical error bars. 
These plots are directly analogous to those shown in Figure 6. 

the method consists of a dynamic measure of the local tur- 
bulent energy on sub-grid scales and sets the flame speed to 



ls^ + C,q 2 , where q is a velocity that characterizes the 

sub-grid turbulence energy content and C, is a constant taken 
to be 4/3. This addition should boost the burning rate during 
the deflagration phase, which will change the density profile 
of the star when the detonation occurs and thus change the 
final yield. It may be our omission of turbulence and flame- 
turbulence interaction models that accounts for the systematic 
difference we see in the yield of IGEs compared to Seitenzahl 
et al. (201 1). While the difference in treatment of the turbulent 
flame is a very important one, it is difficult to decouple from 




Fig. 14. — : Plot of Ms6 Ni vs. age comparing the scaled results 
of this study to the 56 Ni masses inferred from the observations 
of Neill et al. (2009). In red are the points from this study 
with no shift (i.e., tms = Gyr), along with the standard error 
of the mean and a best-fit trend line following the form of 
Equation (5). The green shaded region shows our best-fit line 
with tms = 0.05 - 1 .0 Gyr. In blue are the binned and averaged 
points from Figure 5 of Neill et al. (2009), along with their 
best-fit trend line. 

the dimensionality of the simulations performed. The dimen- 
sionality appears to influence the expansion rate, which may 
or may not be related to the turbulent flame model. Future 
three-dimensional simulations with different turbulent flame 
models will be able to address these uncertainties. 

Seitenzahl et al. (201 1) also mention the choice of />ddt as a 
possible source of the difference between their study and our 
earlier results. As they describe, Seitenzahl et al. (2011) use 
a dynamic measure that calculates the probability of a DDT 
based on the turbulent intensity, density, and fraction of fuel. 
As they note, this prescription is significantly different from 
just fixing pddt and results in fewer DDTs. Considering Ta- 
ble 1 of Seitenzahl et al. (201 1), we can see that the density at 
which the DDT occurred appears to be a function of central 
density due to their DDT criteria: other than their three mod- 
els with only five DDT points, their effective DDT density 
increases with central density. A higher DDT density implies 
an earlier detonation, which means that the star has expanded 
less prior to the detonation, allowing the detonation to burn a 
larger fraction of the star. So it may be possible that the trend 
of increasing IGEs seen in Seitenzahl et al. (201 1) comes from 
the change in their effective DDT density, instead of directly 
from the changing p C Q. Regardless of whether this difference 
in methodology fully explains the difference in results, we can 
state that the differences in methodologies make direct com- 
parison between the two results difficult. 

Finally, the relatively small number of simulations per- 
formed by Seitenzahl et al. (201 1) may not allow them to de- 
termine statistically meaningful trends. Unfortunately, the ex- 
pense of three-dimensional simulations make extensive stud- 
ies difficult. The twelve simulations they performed are a re- 
markable achievement, but our results suggest that it is nec- 
essary to have a distribution of ignition conditions that repro- 
duces the distribution of observed yields. Our statistical en- 
semble does this in two dimensions, but demonstrating such 
a distribution in three dimensions would require significant 
computational resources. 

As Seitenzahl et al. (2011) mention, Meakin et al. (2009) 
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explored the effect of varying the offset from the center of 
the initially burned region (what we call the "match head") 
in the GCD scenario. Their Figure 12 plots mass of NSE el- 
ements (which we have been referring to as IGEs), mass at 
p > 10 7 g cm" 3 , and mass of 56 Ni as functions of central den- 
sity at the time of the detonation; all of their models have the 
same central density at the ignition of the deflagration. We 
note that the GCD scenario is significantly different from the 
DDT scenario that we and Seitenzahl et al. (2011) investi- 
gated, with the principal difference arising in the expansion 
of the star during the deflagration phase of the evolution. In 
the GCD scenario, the star expands during the deflagration as 
the bubbles rise because the displaced mass softens the gravi- 
tational potential. The effect is that GCD models typically ex- 
perience considerably less expansion than DDT models. Also, 
unlike central ignitions, there is very little material burned 
until the detonation occurs. A direct comparison, therefore, 
between the results of Meakin et al. (2009), with yields pre- 
sented as functions of the central density at the ignition of 
the detonation, and our work, considering yields as a function 
of central density at the initial ignition of the deflagration, is 
at best uncertain, especially given that our results show that 
the effect of neutronization is greatest during the deflagration 
phase as illustrated in Figure 7. 

Ropke et al. (2006) performed a set of simulations of SNela 
assuming the pure deflagration model. They found that as 
the central density of the progenitor is increased, the mass of 
56 Ni produced also increases. Their discussion of why they 
produce more iron-group material also applies to our results. 
The yields during the deflagration and detonation phases were 
presented in Figure 7. We find the yield of IGEs increases 
with central density during the deflagration phase. We do not, 
however, find the same increase in 56 Ni mass due to increased 
neutronization. On this point our results disagree with the 
pure deflagration model as studied by Ropke et al. (2006). 
Bravo et al. (1990), however, studied the effects of neutron- 
ization and found that at higher central densities the yield of 
56 Ni decreases, in agreement with our results. 

Similarly to the Seitenzahl et al. (201 1) work, Ropke et al. 
(2006) and Bravo et al. (1990) only performed a modest num- 
ber of simulations. In light of our findings in Section 4.2, we 
note that it may be difficult to draw statistically significant re- 
sults from a single realization. Ropke et al. (2006) and Bravo 
et al. (1990) both appear to use only a single morphology for 
each parameter set, which our results suggest could result in 
deriving an incorrect trend. 

6. SUMMARY AND CONCLUSIONS 

This paper builds on the results presented in Krueger et al. 
(2010), giving more detail of the study and extending the anal- 
ysis. In that paper we showed that, in our 2-d simulations, a 
higher central density of the progenitor star does not impact 
the production of IGEs, but leads to greater neutronization, 
resulting in the production of less 56 Ni. We also discussed the 
relation between the age of the progenitor and the central den- 
sity (see, e.g., Lesaffre et al. 2006), and the relation between 
the brightness and the mass of 56 Ni produced. Thus the state- 
ment that a higher density leads to less 56 Ni is equivalent to 
the statement that an older progenitor will produce a dimmer 
SNIa. In this work we expand on the discussion of Krueger 
et al. (2010) to give more detail of our models and to improve 
upon the age-brightness relation predicted by our simulations. 
In particular, we show that by adding a main-sequence life- 
time to the cooling time our brightness-age relation is steeper, 



more closely matching the observed behavior of older SNela. 

In comparing with other theoretical work, we see that the 
variation of 56 Ni mass with progenitor central density is not a 
settled question. In this paper we further developed the idea 
that, due to the strong nonlinearities of the processes in SNela, 
a statistical study of an ensemble of SNIa simulations may be 
necessary to determine the true trends. For our simulations, 
we find that 15 realizations (morphologies of the initial flame 
surface) are sufficient to characterize the mean trends from 
our models. 

We find that the inner region (out to an enclosed mass of 
0.8 - 1.0 Mq) of the remnant is dominated by 56 Ni. However, 
the stable (non- 56 Ni) IGEs tend to be in "clumps", instead of 
well-mixed throughout this region. This may give rise to 56 Ni 
holes with little or no 56 Ni, depending on the line of sight 
through a SNIa remnant. The outer region of the remnant will 
have more intermediate- and low-mass elements, as the burn- 
ing becomes less efficient for lower densities. As the central 
density increases, the mean 56 Ni mass fraction in the inner re- 
gion drops (roughly 0.8 for p c ,o = 1 x 10 9 g cm" 3 to roughly 
0.6 for p c x) = 5 x 10 9 g cm" 3 ). However, the extent of this re- 
gion (in enclosed-mass space) does not significantly change. 
Variations in the central density affect the sharpness of the 
edges of the stable IGE clumps: a higher central density leads 
to clumps of stable elements that are more sharply defined, as 
well as less mixing between the 56 Ni and the stable IGEs. 

To better connect to observations, we discussed how to dis- 
tinguish the relative ages of SNela with the same brightness 
(in other words, the relative initial central densities of SNela 
that produce the same mass of 56 Ni). We found that, in our 
models, the best measure of the central density is the mass of 
stable IGEs, where higher central density progenitors produce 
more stable IGEs due to their greater rate of neutronization 
during the subsonic deflagration phase. 

We found that a higher central density leads to a shorter 
deflagration phase. Since the rate of neutronization is signif- 
icantly boosted, the total neutronization is greater at a higher 
central density despite there being less time to neutronize. 
The time between the ignition of the first detonation and the 
cessation of burning is independent of central density. 

As noted in Krueger et al. (2010) and described above, our 
choice for the DDT transition density led to an overproduc- 
tion of 56 Ni. In Section 5.1 we provided a recalibration of 
this overall brightness normalization to extrapolate our results 
to an expected average brightness. Future work will be im- 
proved by a better choice of DDT density and we will report 
any quantitative changes to the trends reported here. Due to 
the fundamentally 3-d nature of some of the phenomena in 
a SNIa (such as the turbulent velocities), we plan to extend 
this work by performing 3-d simulations. Because 3-d sim- 
ulations are much more computationally expensive than the 
corresponding 2-d simulations, a study with 3-d simulations 
will by necessity be constrained to a smaller number of simu- 
lations. The choices will be motivated by this work and seek 
to span the parameter space explored here. 
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APPENDIX 

A. RANDOMIZED INITIAL CONDITIONS 
The initial flame surfaces for each realization are determined by spherical harmonics. The formula used is: 

r=150km+J2 A eY"\ (Al) 

i 

where Y" 1 are the spherical harmonics and the coefficients Ag are given in Table Al. A reference implementation is available 
online at http://astronomy.ua.edu/townsley/code/. This suite of simulations used an initial seed value of 1866936915. 

B. RECALIBRATION OF 56 Ni YIELDS 

As discussed above near the beginning of Section 4, on average our suite of simulations exhibited an overproduction of IGEs 
and 56 Ni. The yields of IGEs and 56 Ni are related to the expansion of the WD prior to the ignition of the detonation: as the star 
expands during the deflagration phase, the density decreases and, therefore, less fuel (C and O) remains at sufficiently high density 
to burn to IGEs during the detonation. Thus, the amount of expansion of the WD prior to the DDT principally determines the yield 
of IGEs and 56 Ni (Nomoto et al. 1984; Khokhlov 1991). The amount of expansion is determined by DDT conditions (Jackson 
et al. 2010) along with the energy deposition history, which depends on initial conditions, dimensionality, and the burning model. 

Given our models and the simulations we have performed, the most straightforward way to systematically correct the over- 
production of IGEs is to modify the DDT conditions and repeat the simulations. The value of pddt should be determined by 
the physics of the DDT, but, as the physics of DDTs is incompletely understood and most likely could not be resolved in these 
simulations, we treat pddt as a parameter. The DDT transition is implemented by burning small regions ahead of rising plumes, 
when the tip of the plume reaches pddt- As demonstrated in Jackson et al. (2010), lowering the value of Pddt is tantamount to 
increasing ?ddt because the material that undergoes a DDT is at a higher radius and the time required for a rising bubble of burn- 
ing material to reach that radius is longer. This longer ?ddt results in more expansion of the WD prior to detonation and, because 
the amount of material at higher density decreases, the detonation produces a lower yield of IGEs. In order to address the impact 
of Pddt on our results, the ideal solution would be to repeat our simulations with pddt set to a lower value and therefore obtain 
a longer foDT, more expansion, and more realistic yields. Using the data available from this study and the study of Jackson et al. 
(2010), we can estimate the results of performing a new suite of simulations using a different value of pddt and use these findings 
to inform future studies. We consider two issues here. First we attempt to extrapolate to lower pddt based on the expansion 
characteristics of our models to confirm that for a single value of pddt, the Mice appears to continue to be independent of p C Q. 
Second, we use this justified assumption to extrapolate our yields based on a delay of the DDT transition. 

First we consider how we expect Mice might change if a lower pddt were assumed. Townsley et al. (2009) demonstrates a 
correlation between Mice at ?ige and the mass of all material with a density greater than 2 x 10 7 g cm" 3 at time foDT (which we 
represent by M^T,)- Therefore we can choose a new value of pddt and use this correlation to estimate the final mass of IGEs in 
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TABLE Al: Amplitudes of spherical harmonic perturbations by realization. 



amplitude of perturbation (km) 
Ai 2 A 13 A M Ai 5 A 16 



1 


-50.88 


-51.22 


22.23 


9.979 


-56.23 


2 


-31.03 


51.57 


-43.21 


-55.76 


-30.22 


3 


-39.89 


-28.06 


31.48 


6.848 


-29.80 


4 


-48.89 


-44.21 


1.509 


52.44 


29.00 


5 


29.67 


-74.12 


-53.44 


13.83 


-51.24 


6 


-26.95 


5.169 


-45.79 


-46.34 


14.93 


7 


12.67 


-19.76 


-90.22 


-2.312 


-20.83 


8 


-10.23 


1.206 


-43.74 


-25.57 


10.26 


9 


30.25 


69.91 


-29.37 


-74.17 


-29.50 


10 


-15.76 


-4.818 


14.53 


1.979 


-15.23 


11 


-55.91 


47.08 


-15.86 


-34.44 


-43.00 


12 


1.021 


-0.685 


-33.13 


-48.32 


-22.22 


13 


-10.41 


-32.27 


5.236 


27.06 


-12.46 


14 


-73.60 


-38.16 


-22.76 


41.81 


25.76 


15 


-74.20 


-3.179 


-30.20 


-36.56 


-21.16 


16 


-86.84 


16.28 


19.75 


14.43 


4.924 


17 


12.87 


44.39 


-12.17 


-34.34 


-67.18 


18 


-13.88 


-62.62 


-35.93 


-41.57 


-81.32 


19 


-31.09 


-49.91 


2.953 


57.76 


-5.397 


20 


51.43 


-44.30 


-41.00 


26.83 


-42.83 


21 


-28.15 


11.27 


-37.43 


85.20 


-46.84 


22 


53.67 


4.730 


-52.71 


-46.66 


-65.99 


23 


-21.66 


-12.31 


53.21 


14.14 


-34.18 


24 


-16.25 


-68.30 


-43.63 


111.6 


-33.86 


25 


-56.50 


-31.25 


44.51 


0.098 


-40.66 


26 


-98.62 


-31.80 


-19.40 


45.75 


-3.815 


27 


-39.29 


2.525 


-32.95 


27.10 


-14.95 


28 


29.24 


79.73 


-30.80 


-19.18 


-96.98 


29 


-14.08 


46.98 


-89.82 


-11.95 


28.21 


30 


-24.91 


33.12 


-46.75 


-10.40 


-3.940 



order to verify that changing pddt does not change the fact that the mass of IGEs is independent of p c .o. In order to refine the 
relation between Mice and M®^, we make use of data from this study and data from the study presented in Jackson et al. (2010), 
which varied Pddt- This will enable us to calibrate the relationship between Mjge and for differences in p c o and Pddt- 

The relationship between Mice and M° ^T, appears to be approximately linear, so we assume that Mice = m + ^- As can ^ e 

seen in Figure Bl, the relationship between M°°T, and ^ige depends on p c .o and pddt, so the slope and intercept are allowed to 
be functions of p c ,o and pddt- We test polynomials of varying degree and find that the best functions are quadratic in p Ci o and 
Pddt for both the slope and the intercept. The equations for this fitting are 

M 1G e = m(p c . , pddt) M™ 2 + b{p cXh p DDT ) (B la) 

w(Pc,0,PDDT) = mo + llPcQ + ^plfi + SlPDDT + SlPbuT (Bib) 
b(Pcfi,PDDT) = bo + £lPr,0 + £2Pc,0 + ClPDDT + C2PDDT- ( B 1 C ) 

Fitting this 10-parameter function to the data from this study and the study of Jackson et al. (2010) yields the parameters given 
in Table Bl. As a check, we use this relation at the same pddt as was used in our simulations (10 7 1 g cm" 3 ) and compare the 
estimated values of Mice to the actual values. This is shown in Figure B2. 

With this function, we can extract M™j assuming different values of pddt- Based on the data available, we cannot choose a 
value of pddt lower than what used in this study, but we can test higher values of pddt- Figure B3 shows that for three sample 
values of /?ddt, ^ige is still independent of p c xu formally, we say that Mice is independent of p c q when the magnitude of the 
slope of the line relating these two quantities is less than the uncertainty in the slope. We therefore assume that, even for lower 
values of pddt, Mice is independent of p c u . This allows us to apply a uniform shift in MiGE(constant for all simulations) to correct 
for the overproduction of IGEs. Observations tell us that the mass of IGEs should lie approximately in the range of 0.7-0.9 M Q 
(see, e.g., Woosley et al. 2007), so we can choose to force our mean mass of IGEs to be in the center of this range (0.8 M Q ). 

Now that we have established a reasonable expectation that Mice will continue to be independent of p c o even at other Pddt 
values, we proceed to extrapolate our yields based around this premise. We can relate this change in Mice to a change in the mass 
of 56 Ni through adjustments to ?ddt- The value of foDT, which depends on pddt, is not a free parameter that can be independently 
adjusted in our models; however, in order to correct the mass of Ni we can treat ?ddt as a parameter. Figure B4 shows the 
masses of IGEs and 56 Ni plotted as functions of ?ddt, where each point is a single simulation colored by p c rj. Also shown is a 
trend line based on Equation (B2), explained below. Each of these quantities shows five separate trends, one for each value of 
Pc.o- Performing independent linear fits for each p Cj o gives lines that appear to be correlated: the five lines intersect each other 
near ?ddt = 0, and the slopes appear to be a function of p C Q. Based on these correlations, we derived a new fit using the function 



y(Pc,0,fDDT) = y0 + s(p c fi) fDDT, 
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0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 



M P7 >2 at t DDT ( m o> m p 7 >2 at t DDT (M ) 

Fig. B 1 . — : Relationship between and Mice, demonstrating the dependence of this relationship on p c o and pddt- The data 

in the left panel comes from the study of Jackson et al. (2010), which varied pddt with a constant p c p; data points and trend lines 
are colored by pddt- The data in the right panel comes from this study; data points and trend lines are colored by p c $. 

TABLE Bl: Fit parameters for Equation (Bl). These values assume that p c o is in 10 9 g cm" 3 , pddt is in 10 7 g cm -3 , and both 
masses are in M Q . 



parameter 


value 


mo 


9.099 X KT 1 


7i 


1.453 X 10" 1 


72 


-2.429 x KT 2 


5, 


-4.225 x ltr 1 


Si 


8.555 X 1(T 2 


ho 


8.358 X 10 -3 


£1 


-1.847 x 10 _1 


£2 


2.869 x 1(T 2 


Cl 


6.231 x 10-' 


C2 


-1.196 x KT 1 




0.7 - 

0.6 L- 1 1 1 1 1 1 1 

0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 

M IGE< M p 7 >2. Pc,0> PDDT) (M s ) 

Fig. B2. — : Demonstration of the results of Equation (Bl), showing the comparison of the actual Mice data from simulations to 
the values of Mige estimated from this relation (assuming pddt = 10 7 1 g cm" 3 ). Data points are colored by the study they were 
taken from: blue from this study, red from Jackson et al. (2010). The dotted line shows a perfect correspondence between the 
two measures of Mice- 

where y is the mass of either IGEs or 56 Ni. We found that the best form for the slope s(p c o) is: 

s{p c .o) = a p 2 C Q + b p c .o + c. (B3) 

Minimizing the xt f° r this 4-parameter function results in the parameter values shown in Table B2. 

Given the shift in the mass of IGEs, we can then use the Mtge^cO^ddt) relation from Equations (B2) and (B3) to compute a 
corresponding change in ?ddt- However, since the relation (in particular the slope) depends on p c ,o, the time shift will also depend 
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actual (10 7 ' g cm' 3 ) 
estimated (10', g cm',) 
estimated (10'g g cm "p 
estimated (10 g cm' ) 



p c0 (10 9 gcm' 3 ) 



Fig. B3. — : Demonstration of the results of Equation (Bl), showing the trend of Mice with p c Q, using data from this study. 
The values of Mice from the simulations are shown in red, with green, blue, and magenta showing the estimated Mice assuming 
Pddt = 10 7 1 , 10 7 2 , and 10 7 3 g cm" 3 respectively. The data points are slightly offset horizontally for clarity. This shows that Mice 
is independent of p c> o even if Pddt varies. This figure also demonstrates that a lower value of p DDT leads to a lower mean value 

Of MiGE ■ 




time (s) time (s) 

Fig. B4. — : Yield of IGEs (left panel) and 56 Ni (right panel) as functions of Tddt for the 149 simulations. Also shown are the 
best-fit trend lines in the form given by Equations (B2) and (B3). Lines and data points are colored by p c xy. 1.0 x 10 9 g cm" 3 
(red), 2.0 x 10 9 g cm" 3 (green), 3.0 x 10 9 g cm" 3 (blue), 4.0 x 10 9 g cm" 3 (magenta), 5.0 x 10 9 g cm" 3 (cyan). 

TABLE B2: Fit parameters for Equations (B2) and (B3). These values assume that ?ddt is in seconds, p c is in 10 9 g cm" 3 , and 
mass is in M Q . 



parameter 


MiGE 


M 5( , Ni 


yo 


1.813 


1.550 




1.716 x ur 2 


1.422 x 10~ 2 


b 


-1.859 x 10" 1 


-2.007 x 10-' 


c 


-3.180 x 10-' 


-2.114 x 10-' 



on p c f)\ i.e., we will have A?ddt(Pco)- The five values of A?ddt are given in Table B3. Using this density-dependent time shift 
and the M56 Ni (p Ci o,fDDT) relation, we can then derive a consistent adjustment to the mass of 56 Ni. This adjusted mass of 56 Ni is 
used in Section 5.1. 

C. SIMULATION RESULTS 

The structure and composition of the progenitor models are discussed in detail in Section 2.1, along with the creation of a 
statistical ensemble. The simulations and the code are discussed in Section 2. As we have stressed, the principal result from 
our simulations is a mass of 56 Ni. Our models reproduce the qualitative description given in Woosley et al. (2007) that produces 
a type la supernova light curve. The ideal solution for comparing our results to observed light would be for us to post-process 
our results with a radiation transfer method to actually calculate light curves. In the absence of that technology, we may infer 
observed properties of supernova events suggested by our results. We report some of those here, including stretch, the scaling of 
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TABLE B3: Time shifts for the recalibration of the 56 Ni and IGE masses. 



p Cl (gem 3 ) 


Ar DDT (s) 


1 X 10 9 


0.4850 


2 X 10 9 


0.3800 


3 x 10 9 


0.3273 


4 x 10 9 


0.2999 


5 x 10 9 


0.2884 



brightness with rapidity of decline that is reported in observational results (Howell et al. 2009; Jha et al. 2007). The tabulated 
data are 

1. pcxy. central density of pre-supernova white dwarf (g cm" 3 ); 

2. r. realization number 10 ; 

3. ?ddt: time of the first deflagration-to-detonation transition (s); 

4. MiceI/ddt): mass of Fe-group elements synthesized by time Tddt (M q ); 

5. M56 Ni (fDDT): mass of 56 Ni synthesized by time Tddt (Mq); 

6. ?ige: time that production of Fe-group elements ceases (s); 

7. Mige^ige): mass of Fe-group elements synthesized by time Jige (M Q ); 

8. M56 Ni (fiGE): mass of 56 Ni synthesized by time ?ige (Mq); 

9. fj$ DX : recalibrated 11 value of Tddt (s); 

10. Mf CE (t? GE ): recalibrated 11 value ofM IGE (fiGE) (M ); 

11. M* 6m (t{ GE ): recalibrated 11 value of M 56Ni (f IGE ) (M ); 

12. s: recalibrated 11 value of stretch 12 ; 

13. r cool : cooling time 12 (yr). 

The initial composition for the progenitor is shown in Table 2. The total mass of the star and the total mass of the initial convective, 
isentropic core are both shown for each progenitor in Table 1. Both of these tables are shown in Section 2.1. 

TABLE C 1 : Data extracted from simulations. 



Pc.O 


r 


'DDT 


Mige('ddt) 


M 5 6 Ni (?DDT) 


flGE 


MigeOige) 


Afi6 Ni (flGE) 


? DDT 


^ige^ige) 


^Ni^GE) 


.V 




(g cirT 3 ) 




(s) 


(M ) 


(M ) 


(s) 


(M Q ) 


(M ) 


(s) 


(M ) 


(M ) 




(yr) 


le+9 


1 


1.439 


1.256e-l 


7.744e-2 


2.046 


1.143e+0 


1.01 le+0 


1.924 


9.067e-l 


8.177e-l 


1.203 


b 


le+9 


2 


1.608 


1.756e-l 


1.210e-l 


2.040 


1.026e+0 


9.042e-l 


2.093 


7.900e-l 


7.112e-l 


1.138 


b 


le+9 


3 


1.809 


1.773e-l 


1.229e-l 


2.203 


9.113e-l 


7.991e-l 


2.294 


6.752e-l 


6.061e-l 


1.066 


b 


le+9 


4 


1.696 


1.766e-l 


1.256e-l 


2.083 


1.014e+0 


8.972e-l 


2.181 


7.783e-l 


7.042e-l 


1.134 


b 


le+9 


5 


1.504 


1.490e-l 


9.677e-2 


2.039 


1.041e+0 


9.166e-l 


1.989 


8.052e-l 


7.237e-l 


1.146 


b 


le+9 


6 


1.456 


1.176e-l 


7.158e-2 


1.959 


1.146e+0 


1.015e+0 


1.941 


9.102e-l 


8.217e-l 


1.205 


b 


le+9 


7 


1.474 


1.313e-l 


8.181e-2 


2.008 


1.094e+0 


9.659e-l 


1.959 


8.575e-l 


7.729e-l 


1.177 


b 


le+9 


8 


1.477 


1.164e-l 


6.855e-2 


2.400 


8.780e-l 


7.687e-l 


a 


a 


a 


a 


a 


le+9 


9 


1.682 


1.637e-l 


1.110e-l 


2.095 


9.935e-l 


8.751e-l 


2.167 


7.575e-l 


6.822e-l 


1.119 


b 


le+9 


10 


2.017 


2.081e-l 


1.495e-l 


2.430 


8.959e-l 


7.837e-l 


2.502 


6.599e-l 


5.907e-l 


1.054 


b 


le+9 


1 1 


1.457 


1.310e-l 


8.43 le-2 


1.878 


1.138e+0 


1.008e+0 


1.942 


9.016e-l 


8.155e-l 


1.202 


b 


le+9 


12 


1.556 


1.099e-l 


6.654e-2 


2.025 


1.146e+0 


1.016e+0 


2.041 


9.098e-l 


8.227e-l 


1.206 


b 



Continued on Next Page. . . 



See Appendix A and Section 2 for more details. 

See Appendix B and Section 5.1 for details of the recalibration. 

See Section 5.2 for details of the derivation of s and t coo i. 
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TABLE Cl - Continued 



Pcfl 


r 


'DDT 


Mige('ddt) 


*fe Ni («DDT) 


flGE 


M IG E(flGE) 


M56 Ni (t IGE ) 


'ddt 


M IGE<>iGE) 


M « Ni ( f I*GE) 


.V 




(g cirT 3 ) 




(s) 


(Mq) 


(M ) 


(s) 


(Mq) 


(Mq) 


(s) 


(Mq) 


(Mq) 




(yr) 


le+9 


13 


2.048 


2.017e-l 


1.429e-l 


2.473 


7.969e-l 


6.920e-l 


2.533 


5.608e-l 


4.990e-l 


0.980 




le+9 


14 


1.600 


1.616e-l 


1.135e-l 


2.015 


1.061e+0 


9.410e-l 


2.085 


8.246e-l 


7.480e-l 


1.161 




le+9 


15 


1.441 


1.169e-l 


7.452e-2 


1.863 


1.121e+0 


9.955e-l 


1.927 


8.850e-l 


8.026e-l 


1.194 




le+9 


16 


1.629 


1.756e-l 


1.198e-l 


2.106 


1.028e+0 


9.049e-l 


2.114 


7.923e-l 


7.120e-l 


1.139 




le+9 


17 


1.721 


1.701e-l 


1.204e-l 


2.136 


1.012e+0 


8.951e-l 


2.206 


7.757e-l 


7.021e-l 


1.132 




le+9 


18 


1.407 


1.160e-l 


7.173e-2 


1.906 


1.154e+0 


1.022e+0 


1.892 


9.175e-l 


8.295e-l 


1.210 




le+9 


19 


1.691 


1.690e-l 


1.157e-l 


2.068 


9.690e-l 


8.524e-l 


2.176 


7.330e-l 


6.594e-l 


1.103 




le+9 


20 


1.643 


1.910e-l 


1.341e-l 


2.166 


9.528e-l 


8.358e-l 


2.128 


7.168e-l 


6.428e-l 


1.092 




le+9 


21 


1.299 


9.447e-2 


5.357e-2 


1.824 


1.187e+0 


1.052e+0 


1.784 


9.513e-l 


8.589e-l 


1.226 




le+9 


22 


1.535 


1.459e-l 


9.739e-2 


2.032 


9.850e-l 


8.697e-l 


2.020 


7.489e-l 


6.767e-l 


1.115 




le+9 


23 


1.759 


1.648e-l 


1.131e-l 


2.234 


9.643e-l 


8.498e-l 


2.244 


7.283e-l 


6.568e-l 


1.102 




le+9 


24 


1.666 


1.851e-l 


1.265e-l 


2.179 


8.789e-l 


7.662e-l 


2.151 


6.428e-l 


5.732e-l 


1.041 


k 


le+9 


25 


1.789 


1.865e-l 


1.315e-l 


2.243 


8.950e-l 


7.846e-l 


2.274 


6.590e-l 


5.916e-l 


1.055 


b 


le+9 


26 


1.446 


1.348e-l 


9.107e-2 


1.844 


1.112e+0 


9.889e-l 


1.931 


8.758e-l 


7.959e-l 


1.190 


b 


le+9 


27 


1.392 


9.724e-2 


6.110e-2 


1.845 


1.175e+0 


1.046e+0 


1.877 


9.387e-l 


8.530e-l 


1.223 


b 


le+9 


28 


1.342 


1.132e-l 


7.015e-2 


1.828 


1.096e+0 


9.702e-l 


1.827 


8.598e-l 


7.772e-l 


1.179 


b 


le+9 


29 


1.624 


1.738e-l 


1.169e-l 


2.176 


1.014e+0 


8.911e-l 


2.109 


7.783e-l 


6.981e-l 


1.130 


b 


le+9 


30 


1.571 


1.277e-l 


8.167e-2 


2.000 


1.119e+0 


9.920e-l 


2.056 


8.828e-l 


7.990e-l 


1.192 


b 


2e+9 


1 


1.274 


2.186e-l 


1.056e-l 


1.772 


9.744e-l 


7.991e-l 


1.654 


7.384e-l 


5.879e-l 


1.052 


7.908e+7 


2e+9 


2 


1.338 


2.363e-l 


1.183e-l 


1.827 


8.956e-l 


7.246e-l 


1.718 


6.596e-l 


5.133e-l 


0.992 


7.908e+7 


2e+9 


3 


1.373 


2.242e-l 


1.122e-l 


1.867 


9.631e-l 


7.917e-l 


1.753 


7.270e-l 


5.804e-l 


1.046 


7.908e+7 


2e+9 


4 


1.231 


2.044e-l 


9.844e-2 


1.722 


1.120e+0 


9.398e-l 


1.61 1 


8.838e-l 


7.285e-l 


1.149 


7.908e+7 


2e+9 


5 


1.164 


2.030e-l 


9.237e-2 


1.672 


1.053e+0 


8.705e-l 


1.544 


8.166e-l 


6.592e-l 


1.103 


7.908e+7 


2e+9 


6 


1.278 


1.909e-l 


8.627e-2 


1.779 


1.054e+0 


8.780e-l 


1.658 


8.180e-l 


6.668e-l 


1.109 


7.908e+7 


2e+9 


7 


1.237 


1.942e-l 


8.247e-2 


1.628 


1.090e+0 


9.049e-l 


1.617 


8.541e-l 


6.936e-l 


1.127 


7.908e+7 


2e+9 


g 


1.196 


1.727e-l 


6.563e-2 


1.789 


1.094e+0 


9.092e-l 


1.576 


8.577e-l 


6.979e-l 


1.129 


7.908e+7 


2e+9 


9 


1.281 


2.196e-l 


1.018e-l 


1.742 


1.004e+0 


8.231e-l 


1 .66 1 


7.678e-l 


6.118e-l 


1 .070 


7.908e+7 


2e+9 


10 


1.600 


2.540e-l 


1.408e-l 


2.141 


8.565e-l 


6.948e-l 


1.980 


6.205e-l 


4.835e-l 


0.967 


7.908e+7 


2e+9 


1 1 


1.195 


1.883e-l 


8.399e-2 


1.598 


1.096e+0 


9.183e-l 


1.575 


8.598e-l 


7.070e-l 


1.135 


7.908e+7 


2e+9 


12 


1.134 


1.356e-l 


4.169e-2 


1.604 


1.194e+0 


1.007e+0 


1.514 


9.575e-l 


7.957e-l 


1.190 


7.908e+7 


2e+9 


13 


1.524 


2.363e-l 


1.168e-l 


1.957 


8.689e-l 


6.994e-l 


1.904 


6.328e-l 


4.882e-l 


0.971 


7.908e+7 


2e+9 


14 


1.242 


2.062e-l 


1.029e-l 


1.727 


1.016e+0 


8.468e-l 


1.622 


7.795e-l 


6.356e-l 


1.087 


7.908e+7 


2e+9 


15 


1.024 


1.429e-l 


5.202e-2 


1.569 


1.170e+0 


9.880e-l 


1.404 


9.340e-l 


7.768e-l 


1.179 


7.908e+7 


2e+9 


16 


1.347 


2.412e-l 


1.257e-l 


1.816 


9.349e-l 


7.629e-l 


1.727 


6.988e-l 


5.516e-l 


1.024 


7.908e+7 


2e+9 


17 


1.275 


2.046e-l 


9.945e-2 


1.704 


1.058e+0 


8.838e-l 


1.655 


8.223e-l 


6.726e-l 


1.112 


7.908e+7 


2e+9 


18 


1.019 


1.490e-l 


5.379e-2 


1.586 


1.189e+0 


1.003e+0 


1.399 


9.533e-l 


7.921e-l 


1.188 


7.908e+7 


2e+9 


19 


1.206 


1.913e-l 


8.005e-2 


1.619 


1.105e+0 


9.195e-l 


1.586 


8.685e-l 


7.083e-l 


1.136 


7.908e+7 


2e+9 


20 


1.298 


2.323e-l 


1.153e-l 


1.716 


1.013e+0 


8.333e-l 


1.678 


7.774e-l 


6.220e-l 


1.077 


7.908e+7 


2e+9 


21 


0.994 


1.383e-l 


4.546e-2 


1.503 


1.225e+0 


1.034e+0 


1.374 


9.892e-l 


8.231e-l 


1.206 


7.908e+7 


2e+9 


22 


1.197 


1.861e-l 


8.048e-2 


1.792 


1.043e+0 


8.660e-l 


1.577 


8.065e-l 


6.548e-l 


1.100 


7.908e+7 


2e+9 


23 


1.337 


2.027e-l 


9.030e-2 


1.780 


9.984e-l 


8.230e-l 


1.717 


7.624e-l 


6.117e-l 


1.070 


7.908e+7 


2e+9 


24 


1.237 


2.299e-l 


1.069e-l 


1.870 


8.652e-l 


6.905e-l 


1.617 


6.291e-l 


4.792e-l 


0.963 


7.908e+7 


2e+9 


25 


1.338 


2.190e-l 


1.059e-l 


1.736 


9.844e-l 


8.099e-l 


1.718 


7.484e-l 


5.986e-l 


1.060 


7.908e+7 


2e+9 


26 


1.135 


1.792e-l 


7.881e-2 


1.557 


1.107e+0 


9.315e-l 


1.515 


8.714e-l 


7.202e-l 


1.144 


7.908e+7 


ze+y 


97 


1 119 
1 . 1 1Z 


1 99Qo 1 


A Q 1 8p> 9 

t.yioe-z 


1 S^Q 


1 181 t*j-C\ 


1 .UUoc+U 


1 AQ9 




7 Q7(V 1 
/ .y /Uc-1 


1 1 Q1 


7 QflS^-u7 


2e+9 


28 


1.014 


1.529e-l 


5.527e-2 


1.451 


1.142e+0 


9.583e-l 


1.394 


9.059e-l 


7.471e-l 


1.161 


7.908e+7 


2e+9 


29 


1.193 


2.064e-l 


8.958e-2 


1.766 


1.021e+0 


8.368e-l 


1.573 


7.845e-l 


6.255e-l 


1.080 


7.908e+7 


2e+9 


30 


1.081 


1.509e-l 


5.572e-2 


1.529 


1.196e+0 


1.008e+0 


1.461 


9.596e-l 


7.971e-l 


1.191 


7.908e+7 


3e+9 


1 


1.111 


2.392e-l 


8.266e-2 


1.533 


1.026e+0 


7.969e-l 


1.438 


7.899e-l 


5.726e-l 


1.040 


2.879e+8 


3e+9 


2 


1.064 


2.459e-l 


8.436e-2 


1.503 


1.058e+0 


8.240e-l 


1.391 


8.217e-l 


5.997e-l 


1.061 


2.879e+8 


3e+9 


3 


1.285 


2.507e-l 


8.898e-2 


1.723 


8.638e-l 


6.488e-l 


1.612 


6.278e-l 


4.244e-l 


0.911 


2.879e+8 



Continued on Next Page. 
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TABLE Cl - Continued 



Pc.o r 


'DDT 


Mige('ddt) M 5 6 Ni («DDT) 


flGE 


MlQE(flGE) M 5 6 Ni («IGE) 


'ddt 


M I*GE( f lGE) M « Nl ( f I*GE) s 




(g cirT 3 ) 


(s) 


(M Q ) (M ) 


(s) 


(M Q ) (Mq) 


(s) 


(M ) (M ) 


(yr) 



3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
3e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 
4e+9 



10 
1 I 

12 
13 
14 
15 
16 
17 
IS 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
1 1 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 



1.174 
1.022 
1.024 
0.998 
1.039 
1.141 
1.401 
1.046 
0.965 
1.341 
1.076 
0.947 
1.129 
1.129 
0.890 
1.204 
1.176 
0.895 
1.028 
1.213 
1.038 
1.203 
0.984 
1.003 
0.921 
1.011 
1.038 
0.996 
0.958 
1.182 
1.092 
0.898 
0.935 
0.943 
1.041 
1.009 
1.306 
0.966 
0.911 
1.243 
0.991 
0.857 
1.072 
0.982 
0.761 
1.103 
1.087 
0.780 
0.954 
1.130 
0.941 
1.073 
0.922 



2.375e-l 
2.249e-l 
1.950e-l 
2.130e-l 
2.016e-l 
2.510e-l 
2.650e-l 
2.204e-l 
1.622e-l 
2.612e-l 
2.186e-l 
1.748e-l 
2.631e-l 
2.272e-l 
1.738e-l 
2.529e-l 
2.544e-l 
1.713e-l 
2.095e-l 
2.392e-l 
2.525e-l 
2.463e-l 
2.020e-l 
1.638e-l 
1.909e-l 
2.293e-l 
1.857e-l 
2.532e-l 
2.501e-l 
2.528e-l 
2.416e-l 
2.296e-l 
2.141e-l 
2.383e-l 
2.287e-l 
2.551e-l 
2.604e-l 
2.429e-l 
1.863e-l 
2.776e-l 
2.307e-l 
1.973e-l 
2.842e-l 
2.230e-l 
1.836e-l 
2.649e-l 
2.577e-l 
1.877e-l 
2.283e-l 
2.532e-l 
2.751e-l 
2.644e-l 
2.217e-l 



8.843e-2 
7.25 le-2 
4.817e-2 
5.642e-2 
5.049e-2 
8.595e-2 
1.099e-l 
6.880e-2 
3.069e-2 
9.459e-2 
7.401e-2 
4.569e-2 
9.851e-2 
7.756e-2 
4.342e-2 
8.665e-2 
9.120e-2 
4.005e-2 
6.057e-2 
7.568e-2 
8.053e-2 
9.050e-2 
5.757e-2 
3.844e-2 
5.010e-2 
6.818e-2 
4.406e-2 
6.652e-2 
5.462e-2 
5.884e-2 
6.524e-2 
4.362e-2 
3.618e-2 
4.479e-2 
3.968e-2 
5.758e-2 
6.975e-2 
5.893e-2 
2.445e-2 
7.328e-2 
5.550e-2 
3.404e-2 
8.135e-2 
5.106e-2 
3.043e-2 
6.228e-2 
6.216e-2 
2.909e-2 
4.964e-2 
5.409e-2 
5.869e-2 
6.909e-2 
4.436e-2 



1 .75 1 

1.556 
1.531 
1.414 
1.538 
1.652 
1.881 
1.484 
1.478 
1.777 
1.578 
1.413 
1.664 
1.589 
1.426 
1.639 
1.608 
1.403 
1.509 
1.626 
1.647 
1.629 
1.385 
1.433 
1.375 
1.490 
1.473 
1.449 
1.376 
1.581 
1.683 
1.344 
1.436 
1.390 
1.524 
1.513 
1.739 
1.385 
1.411 
1.693 
1.432 
1.295 
1.520 
1.475 
1.272 
1.518 
1.549 
1.286 
1.456 
1.637 
1.506 
1.538 
1.397 



9.808e-l 
1.032e+0 
1.110e+0 
1.120e+0 
1.104e+0 
9.048e-l 
8.474e-l 
1.054e+0 
1.216e+0 
7.783e-l 
1.005e+0 
1.159e+0 
8.836e-l 
1.020e+0 
1.180e+0 
9.155e-l 
9.708e-l 
1.195e+0 
1.058e+0 
9.500e-l 
8.788e-l 
9.745e-l 
1.109e+0 
1.172e+0 
1.123e+0 
1.025e+0 
1.137e+0 
1.050e+0 
1.083e+0 
9.215e-l 
9.218e-l 
1.095e+0 
1.109e+0 
l.llle+0 
1.029e+0 
9.731e-l 
8.840e-l 
1.078e+0 
1.202e+0 
8.450e-l 
1.040e+0 
1.174e+0 
9.316e-l 
1.068e+0 
1.224e+0 
9.185e-l 
9.537e-l 
1.228e+0 
1.019e+0 
8.991e-l 
8.969e-l 
9.503e-l 
1.097e+0 



7.627e-l 
8.054e-l 
8.787e-l 
8.818e-l 
8.693e-l 
6.815e-l 
6.409e-l 
8.297e-l 
9.813e-l 
5.657e-l 
7.881e-l 
9.341e-l 
6.622e-l 
8.019e-l 
9.478e-l 
6.899e-l 
7.453e-l 
9.655e-l 
8.335e-l 
7.260e-l 
6.479e-l 
7.522e-l 
8.828e-l 
9.547e-l 
8.973e-l 
7.908e-l 
9.113e-l 
7.802e-l 
8.002e-l 
6.633e-l 
6.719e-l 
8.208e-l 
8.331e-l 
8.314e-l 
7.590e-l 
7.033e-l 
6.309e-l 
8.134e-l 
9.346e-l 
5.828e-l 
7.836e-l 
9.107e-l 
6.651e-l 
8.094e-l 
9.556e-l 
6.502e-l 
6.863e-l 
9.640e-l 
7.587e-l 
6.326e-l 
6.151e-l 
6.829e-l 
8.284e-l 



1.501 
1.349 
1.351 
1.325 
1.367 
1.469 
1.728 
1.373 
1.292 
1.669 
1.403 
1.274 
1.456 
1.456 
1.217 
1.531 
1.503 
1.222 
1.355 
1.540 
1.365 
1.530 
1.312 
1.331 
1.249 
1.339 
1.366 
1.295 
1.258 
1.482 
1.392 
1.198 
1.235 
1.243 
1.341 
1.309 
1.606 
1.266 
1.211 
1.543 
1.291 
1.157 
1.372 
1.282 
1.061 
1.403 
1.387 
1.080 
1.254 
1.430 
1.241 
1.373 
1.222 



7.448e-l 
7.958e-l 
8.739e-l 
8.841e-l 
8.679e-l 
6.687e-l 
6.114e-l 
8.182e-l 
9.799e-l 
5.423e-l 
7.688e-l 
9.226e-l 
6.475e-l 
7.843e-l 
9.438e-l 
6.794e-l 
7.347e-l 
9.594e-l 
8.223e-l 
7.139e-l 
6.428e-l 
7.384e-l 
8.728e-l 
9.360e-l 
8.866e-l 
7.894e-l 
9.013e-l 
8.138e-l 
8.465e-l 
6.855e-l 
6.857e-l 
8.589e-l 
8.732e-l 
8.749e-l 
7.930e-l 
7.370e-l 
6.479e-l 
8.423e-l 
9.661e-l 
6.090e-l 
8.035e-l 
9.383e-l 
6.956e-l 
8.320e-l 
9.875e-l 
6.824e-l 
7.177e-l 
9.919e-l 
7.827e-l 
6.631e-l 
6.608e-l 
7.142e-l 
8.612e-l 



5.383e-l 
5.810e-l 
6.544e-l 
6.575e-l 
6.449e-l 
4.571e-l 
4.165e-l 
6.054e-l 
7.569e-l 
3.414e-l 
5.637e-l 
7.097e-l 
4.378e-l 
5.775e-l 
7.235e-l 
4.655e-l 
5.210e-l 
7.41 le-1 
6.092e-l 
5.017e-l 
4.236e-l 
5.278e-l 
6.584e-l 
7.304e-l 
6.729e-l 
5.665e-l 
6.869e-l 
5.443e-l 
5.642e-l 
4.274e-l 
4.359e-l 
5.848e-l 
5.972e-l 
5.954e-l 
5.231e-l 
4.673e-l 
3.950e-l 
5.775e-l 
6.986e-l 
3.468e-l 
5.477e-l 
6.747e-l 
4.291e-l 
5.735e-l 
7.196e-l 
4.142e-l 
4.504e-l 
7.281e-l 
5.227e-l 
3.966e-l 
3.791e-l 
4.470e-l 
5.924e-l 



1.013 
1.047 
1.100 
1.102 
1.093 
0.943 
0.903 
1.065 
1.167 
0.821 
1.033 
1.137 
0.924 
1.044 
1.146 
0.950 
0.999 
1.157 
1.068 
0.982 
0.910 
1.004 
1.103 
1.150 
1.113 
1.036 
1.122 
1.018 
1.034 
0.914 
0.922 
1.050 
1.059 
1.058 
1.001 
0.952 
0.881 
1.044 
1.130 
0.828 
1.021 
1.114 
0.916 
1.041 
1.143 
0.901 
0.936 
1.149 
1.000 
0.883 
0.864 
0.933 
1.055 



2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
2.879e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
5.682e+8 
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'ddt 


Alien (fmvr ) 
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^DDT 
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{ rr r*m *\ 

Kb Lm / 






(Mr \ 


V ivi ) 




V ivi Q ) 


V ivl Q / 










iyi ) 


4e+9 


27 


0.909 


1.905e-l 


3.106e-2 


1.348 


1.167e+0 


9.091e-l 


1.209 


9.308e-l 


6.731e-l 


1.113 


5.682e+8 


4e+9 


28 


0.873 


2.169e-l 


4.070e-2 


1.321 


1.121e+0 


8.550e-l 


1.173 


8.848e-l 


6.191e-l 


1.075 


5.682e+8 


4e+9 


29 


0.953 


2.442e-l 


4.582e-2 


1.522 


9.938e-l 


7.150e-l 


1.253 


7.577e-l 


4.790e-l 


0.963 


5.682e+8 


4e+9 


30 


0.912 


1.922e-l 


2.870e-2 


1.418 


1.164e+0 


8.932e-l 


1.212 


9.284e-l 


6.573e-l 


1.102 


5.682e+8 


5e+9 


1 


0.993 


2.659e-l 


5.307e-2 


1.441 


9.871e-l 


6.823e-l 


1.282 


7.510e-l 


4.344e-l 


0.921 


9.625e+8 


5e+9 


2 


0.964 


2.683e-l 


4.278e-2 


1.357 


1.040e+0 


7.268e-l 


1.252 


8.037e-l 


4.789e-l 


0.962 


9.625e+8 


5e+9 


3 


1.151 


2.732e-l 


4.979e-2 


1.593 


8.698e-l 


5.775e-l 


1.439 


6.338e-l 


3.296e-l 


0.807 


9.625e+8 


5e+9 


4 


1.010 


2.522e-l 


4.409e-2 


1.589 


9.855e-l 


6.843e-l 


1.298 


7.495e-l 


4.364e-l 


0.923 


9.625e+8 


5e+9 


5 


0.835 


2.397e-l 


3.034e-2 


1.303 


1.073e+0 


7.673e-l 


1.123 


8.366e-l 


5.195e-l 


0.998 


9.625e+8 


5e+9 


6 


0.919 


2.352e-l 


3.124e-2 


1.425 


1.094e+0 


7.91 le-1 


1.207 


8.575e-l 


5.432e-l 


1.017 


9.625e+8 


5e+9 


7 


0.871 


2.459e-l 


3.091e-2 


1.368 


1.071e+0 


7.635e-l 


1.159 


8.346e-l 


5.157e-l 


0.994 


9.625e+8 


5e+9 


8 


1.079 


2.583e-l 


4.057e-2 


1.553 


9.488e-l 


6.520e-l 


1.367 


7.127e-l 


4.042e-l 


0.891 


9.625e+8 


5e+9 


9 


1.028 


2.723e-l 


3.992e-2 


1.438 


9.257e-l 


6.253e-l 


1.316 


6.897e-l 


3.775e-l 


0.862 


9.625e+8 


5e+9 


10 


1.321 


2.654e-l 


3.999e-2 


1.744 


8.197e-l 


5.304e-l 


1.609 


5.836e-l 


2.825e-l 


0.746 


9.625e+8 


5e+9 


1 1 


0.945 


2.598e-l 


4.614e-2 


1.369 


1.020e+0 


7.213e-l 


1.234 


7.844e-l 


4.735e-l 


0.958 


9.625e+8 


5e+9 


12 


0.893 


2.114e-l 


2.626e-2 


1.402 


1.186e+0 


8.934e-l 


1.181 


9.502e-l 


6.455e-l 


1.094 


9.625e+8 


5e+9 


13 


1.168 


2.718e-l 


4.903e-2 


1.658 


8.762e-l 


5.797e-l 


1.457 


6.401e-l 


3.318e-l 


0.810 


9.625e+8 


5e+9 


14 


0.957 


2.477e-l 


4.220e-2 


1.356 


1.051e+0 


7.527e-l 


1.246 


8.148e-l 


5.048e-l 


0.985 


9.625e+8 


5e+9 


15 


0.816 


2.211e-l 


3.014e-2 


1.230 


1.168e+0 


8.686e-l 


1.104 


9.317e-l 


6.207e-l 


1.076 


9.625e+8 


5e+9 


16 


0.962 


2.841e-l 


5.486e-2 


1.476 


9.264e-l 


6.173e-l 


1.250 


6.903e-l 


3.694e-l 


0.854 


9.625e+8 


5e+9 


17 


0.952 


2.412e-l 


3.602e-2 


1.394 


1.017e+0 


7.242e-l 


1.240 


7.806e-l 


4.764e-l 


0.960 


9.625e+8 


5e+9 


18 


0.708 


1.996e-l 


2.283e-2 


1.240 


1.231e+0 


9.288e-l 


0.996 


9.945e-l 


6.809e-l 


1.118 


9.625e+8 


5e+9 


19 


1.056 


2.749e-l 


4.684e-2 


1.514 


9.454e-l 


6.377e-l 


1.345 


7.093e-l 


3.898e-l 


0.876 


9.625e+8 


5e+9 


20 


1.017 


2.628e-l 


3.967e-2 


1.450 


9.610e-l 


6.549e-l 


1.305 


7.250e-l 


4.070e-l 


0.894 


9.625e+8 


5e+9 


21 


0.738 


2.054e-l 


2.201e-2 


1.248 


1.230e+0 


9.259e-l 


1.026 


9.936e-l 


6.781e-l 


1.116 


9.625e+8 


5e+9 


22 


0.872 


2.439e-l 


4.205e-2 


1.420 


1.057e+0 


7.587e-l 


1.160 


8.213e-l 


5.108e-l 


0.990 


9.625e+8 


5e+9 


23 


1.090 


2.600e-l 


3.872e-2 


1.487 


9.235e-l 


6.318e-l 


1.379 


6.874e-l 


3.839e-l 


0.869 


9.625e+8 


5e+9 


24 


0.877 


2.843e-l 


4.293e-2 


1.417 


9.215e-l 


6.055e-l 


1.165 


6.854e-l 


3.577e-l 


0.840 


9.625e+8 


5e+9 


25 


1.034 


2.672e-l 


4.927e-2 


1.469 


9.410e-l 


6.43 le-1 


1.323 


7.050e-l 


3.952e-l 


0.881 


9.625e+8 


5e+9 


26 


0.898 


2.510e-l 


4.255e-2 


1.344 


1.098e+0 


7.913e-l 


1.186 


8.615e-l 


5.435e-l 


1.017 


9.625e+8 


5e+9 


27 


0.985 


2.265e-l 


3.279e-2 


1.413 


1.094e+0 


8.063e-l 


1.274 


8.583e-l 


5.584e-l 


1.029 


9.625e+8 


5e+9 


28 


0.795 


2.264e-l 


2.744e-2 


1.272 


1.146e+0 


8.445e-l 


1.084 


9.096e-l 


5.966e-l 


1.059 


9.625e+8 


5e+9 


29 


0.891 


2.486e-l 


3.129e-2 


1.406 


1.027e+0 


7.174e-l 


1.180 


7.907e-l 


4.695e-l 


0.954 


9.625e+8 


5e+9 


30 


0.868 


2.1626-1 


2.560e-2 


1.343 


1.160e+0 


8.600e-l 


1.156 


9.239e-l 


6.122e-l 


1.070 


9.625e+8 



il As discussed in Section 3.1.2, the results of p c Q = 1.0 X 10' g cm 3 , realization 8 were excluded from the analysis. 

b As discussed in Section 5.2, the results of Lesaffre et al. (2006) do not allow us to generate cooling times for results with p c o = 1 .0 X 10 9 g cm 



